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Abstract. We study the distribution, in the space of Satake parameters, of local components of 
Siegel cusp forms of genus 2 and growing weight k, subject to a specific weighting which allows 
us to apply results concerning Bessel models and a variant of Petersson's formula. We obtain 
for this family a quantitative local equidistribution result, and derive a number of consequences. 
In particular, we show that the computation of the density of low-lying zeros of the spinor L- 
functions (for restricted test functions) gives global evidence for a well-known conjecture of Bocherer 
concerning the arithmetic nature of Fourier coefficients of Siegel cusp forms. 



1. Introduction 

1.1. Motivation. The motivation behind this paper lies in attempts to understand what is a 
correct definition of a family of cusp forms, either on GL(n) or on some other reductive algebraic 
group. The basic philosophy (or strategy) underlying our work is the following form of local-global 
principle: given a "family" II of cusp forms, for any finite place v, the local components % v of the 
elements ir £ II (which are represented as restricted tensor products ir = (g>ir v over all places) should 
be well-behaved, and more specifically, under averaging over finite subsets of the family, (ir v ) should 
become equidistributed with respect to a suitable measure fj, v . Readers already familiar with work 
on families of Dirichlet characters, classical modular forms on GL(2), or families of L-functions 
of abelian varieties, will recognize that this principle is implicit in much of these works, through 
the orthogonality relations for Dirichlet characters, the trace formula (or the Petersson formula) 
and the "vertical" Sato- Tate laws over finite fields. The expected outcome is that, for instance, 
averages over the family of values of L-functions L(so,ir) at some point so, at least on the right 
of the critical line, should be directly related to the Euler product corresponding to local averages 
computed using fj, v . (For a general informal survey of this point of view, see |29j). 

One of our goals is to give an example where this strategy can be implemented in the case of 
holomorphic cusp forms on GSp(4) (i.e., Siegel modular forms), and to derive some applications of 
it. In particular, we will proveo 

Theorem 1.1. For k ^ 2, let §>* k be a Hecke basis of the space of Siegel cusp forms on Sp(4, Z), 
and let § k be the set of those F 6 St which are not Saito-Kurokawa lifts. For F G §* k , let 

F(Z) = ^a(F,T)e(Tr(TZ)) 

T>0 

be its Fourier expansion, where T runs over symmetric positive- definite semi-integral matrices, and 
let 

(i.i.i) <4 = ^r 2k r(k - ink - 2)^f^. 
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Let L(F, s) denote the finite part of the spin L-function of F, an Euler product of degree 4 over 
primes with local factors of the form 

L p (F, s) = (1 - a p p~ s )-\l - b p p- s r\l - a^p-THl ~ V) -1 , %>, \ E (C x ) 2 . 

Then, for any s G C suc/i that Re(s) > 1, we have 

(1.1.2) lim £ ^ fe L(F, S ) = C(s + ±)L( X 4,s + l), 



where ((s) denotes the Riemann zeta function and L(x4,s) is the L-function associated to the 
unique Dirichlet character of conductor 4. For s ^ 3/2, one can replace § 2 k with S^. 

More generally, for all primes p there exist measures p p on (C x ) 2 , which are in fact supported on 
(S 1 ) 2 , with the following property: for any irreducible r-dimensional representation p o/GSp(4, C), 
let L(F, p, s) denote the associated Langlands L-function, an Euler product of degree r ^ 1 over 
primes with local factors of the form 

r 

L p (F,p,s) =Yl(l-Qi(a p ,b p )p- s )-\ op, b p G(C x ) 2 . 

i=l 

where Qi(x,y) is a polynomial in x, y, x^ 1 , y^ 1 . Then for any s E C such that Re(s) > so, with 
so depending on p, we have 

(1.1.3) lim ^ U ^L{F,p, 8 )=]lff[(l-Q i (a,b)p-')- 1 dfi p (a,b), 

FGS* fe V i=l 

where the right-hand side converges absolutely. 

The weight ui^ which is introduced in this theorem is natural because of our main tool, which is 
a (rather sophisticated) extension of the classical Petersson formula to the case of Siegel modular 
forms of genus 2; see Propositions 13.31 and I3T61 In fact, we can work with more general weights 
u kdA ( as defined in the next section) which involve averages of a(F,T) over positive definite T 
with a fixed discriminant. 

The quantitative local equidistribution leads naturally to a result on the distribution of low-lying 
zeros: 

Theorem 1.2 (Low- lying zeros). Let ip : R — » R be an even Schwartz function such that the 
Fourier transform 



0(t) = / <p(x)t 
JR 



has compact support contained in [—a, a], where a < 4/15. For F G S^, assume the Riemann 
Hypothesis: the zeros of L(F,s) in the critical strip < Re(s) < 1 are of the form 

1 

p = - + ^7 

with 7 6 R. Define 

'7 



d <p( f )=J2H 1 log k 



P 



where p ranges over the zeros of L(F,s) on the critical line, counted with multiplicity. 
Then we have 

(1.1.4) lim V ^ k D^{F)= [ <p(x)d* Sp (x), 



where o~sp is the "symplectic symmetry" measure given by 

dasp = dx — ^ , 5q Dime mass at 0. 

This result raises interesting questions concerning the notion of "family" of cusp forms, especially 
from the point of view of the notion of symmetry type that has arisen from the works of Katz- 
Sarnak [25]. Indeed, the limit measure above is the one that arises from symplectic symmetry 
types, i.e., from the distribution of eigenvalues close to 1 of symplectic matrices of large size, when 
renormalized to have averaged spacing equal to 1. In general, it is expected that some cusp forms 
will exhibit this symmetry when some kind of infinite symplectic group occurs as "monodromy 
group" for the family, in the way that generalizes the Chebotarev and Deligne equidistribution 
theorems. 

We do not believe that this is the case here, and rather expect that the limit measure in the 
theorem is due in part to the presence of the weight u^u used in the averages involved. Precisely, 
we expect that the correct symmetry type, without weight, is orthogonal, in the sense that for ip 
with support in ] — 1, 1[, we should havqj 

(1.1.5) j-L D A F ) — > / V>(x)dao(x) 

where dao(x) = dx + ^ is the corresponding measure for eigenvalues close to 1 of orthogonal 
matrices. 

Intuitively, this should be related to the fact that the point 1/2 is a critical special value - 
in the sense of Deligne - for the spin L- functions of cusp forms F £ S£ (whereas 1/2 is not for 
real quadratic characters for example, which are the typical example where symplectic symmetry 
is expected), similar to the special role of the eigenvalue 1 for orthogonal matrices, but not for 
symplectic ones. 

Now the natural question is why should the weight ui^j. have such an effect? (This is especially 
true because it may look, at first, just like an analogue of the weight involving the Petersson norm of 
classical modular forms which has been used very frequently without exhibiting any such behavior, 
e.g., in the works of Iwaniec, Luo and Sarnak [21] and Duenez-Miller jll|.) 

The point is that this weight itself contains arithmetic information related to central L-values 
of the Siegel cusp forms. Indeed, we will see in Section E3] that Theorem II .21 can be interpreted con- 
vincingly - assuming an orthogonal symmetry as in fll . 1 .5 j) - as evidence for a beautiful conjecture 
of Bocherer (see [3] or [191 Intr.]) which suggests in particular a relation of the type 

(1.1.6) \a(F, 1)| 2 ~ L(F, 1/2)L(F x X4 , 1/2) 

(where the ~ sign means equality up to non-zero factors "unrelated to central critical values" ; this 
version of the conjecture is that proposed by Furusawa and Martin [161 §1, (1.4)]). We therefore 
consider that Theorem 11.21 provides suggestive global evidence towards these specific variants of 
Bocherer's conjecture. Note that, at the current time, this conjecture is not rigorously known for 
any cusp form in which is not a Saito-Kurokawa lift. 

Remark 1.3. One can easily present analogues of the phenomenon in Theorem [L2j as we understand 
it, in the setting of random matrices. For instance, if \i n denotes the probability Haar measure on 
S02n(IR)> one may consider the measures 

dv n (g) = c n det(l - g)dfi n (g), 



We do not try to predict whether odd or even orthogonal symmetry should occur; these could be distinguished 
most simply by computing the 2-level density for test functions with restricted support as done in |33| for classical 
modular forms (one can also attempt to study the low- lying zeros for test functions with support larger than ] — 1, 1[, 
as done in [24] . 
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where c n > is the constant that ensures that v n is a probability measure]! The distribution of 
the low- lying eigenvalues of g E S02 n (K), when computed using this measure, will clearly differ 
from that arising from Haar measure (intuitively, by diminishing the influence of matrices with an 
eigenvalue close to 1, the factor det(l — g) will produce a repulsion effect similar to what happens 
for symplectic matrices.) 

Readers familiar with the case of GL(2)-modular forms but not with Siegel modular forms (or 
with their representation-theoretic interpretation) may look at the Appendix where we discuss 
briefly the analogies and significant differences between our results and some more elementary 
GL(2)-versions. 

For orientation, we add the following quick remarks: (1) the spin L-function of F E 2>%, has 
analytic conductor (see [231 P- 95] for the definition) of size k 2 ; (2) the cardinality of (i.e., the 
dimension of S 2 fc) is of order of magnitude k 3 (see, e.g., [271 Cor. p. 123] for the space M 2 fc of all 
Siegel modular forms of weight 2k, and [27J p. 69] for the size of the "correction term" M 2 fc/S 2 fc); 
(3) as already mentioned, the spin L-function is self-dual with functional equation involving the 
sign +1 for all L E S 2fc . 

Apart from the treatment of low-lying zeros, we do not "enter the critical strip" in this paper. 
However, we hope to come back to the problem of extending Theorem 11.11 to averages at points 
inside the critical strip, and we may already remark that, if a statement like (|1.1.6|) is valid, the 
weight already involves some critical values of L-functions (in fact, of an L-function of degree 8. 

1.2. Local equidistribution statement. In order to state our main result on local spectral 
equidistribution of Siegel modular forms, we begin with some preliminary notation concerning 
cuspidal automorphic representations of G(A) = GSp(4, A). 

Let 7r be a cuspidal automorphic representation of G(A), which we assume to be unramified at 
all finite places and with trivial central character. It is isomorphic to a restricted tensor product 
vr = <S>„7r„ where, for all places, ir v is an irreducible admissible unitary representation of the local 
group G(Q„). 

By our assumption ir p is unramified for all primes p and so the natural underlying space for 
local equidistribution at p (when considering families) is the set X p of unramified unitary infinite- 
dimensional irreducible representations of G(Q P ) with trivial central character. This set has a 
natural topology, hence a natural cr-algebra. 

We now proceed quite concretely to give natural coordinates on X p from which the measurable 
structure is obvious. By [7J, any 7r p E X p can be identified with the unique unramified constituent 
of a representation xi x X2 x & induced from a character of the Borel subgroup which is defined as 
follows using unramified (not necessarily unitary!) characters xii X2, & of Q p : 

(a\ * * * \ 



a 2 



Xa^ 1 



Xi(ai)X2(a 2 )<r(A). 



\ * Aa 2 7 

Having trivial central character means that 



X1X2V 2 = 1, 



and since the characters are unramified, it follows that ir p is characterized by the pair (a, b) = 
(cr(p), cr(p)xi(p)) S C* xC*. The classification of local representations of G(Q P ) (see for instance [38l 
Proposition 3.1]) implies that the local parameters satisfy 

(1.2.1) < |CI|, \b\ < yft. 



3 In fact, although this seems irrelevant, a computation with the moments of characteristic polynomials of orthogonal 
matrices shows that c n = 1/2 for all n. 



There are some identifications between the representations associated to different (a, 6), coming 
from the Weyl W group of order 8 generated by the transformations 

(1.2.2) (a, b) i-> (b, a), (a, b) i-> (a -1 , b), (a, 6) i-> (a, b' 1 ). 

We will denote by Y p the quotient of the set of (a, 6) satisfying the upper-bounds (|1.2.ip . modulo 
the action of W. This has the quotient topology and quotient a-algebra, and we identify X p with 
a subset of Y p using the parameters (a, b) described above. We will also denote by X C X p the 
subset of tempered representations; under the identification of X p with a subset of Y p , the set X 
corresponds precisely to \a\ = \b\ = 1. Note that this subset is indeed independent of p. 

In applications to L- functions, the local- global nature of automorphic representations is reflected 
not only in the existence of local components, but in their "independence" (or product structure) 
when p varies. To measure this below, we will also need to consider, for any finite set of primes S, 
the maps 

which have image in the space 

x s = n x p 

pes 

and can be identified with a subset of 

Ys=ilY p . 

pes 

Now we come back to Siegel modular forms. Let §k = §fc(Sp(4, Z)) be the space of Siegel 
cusp forms of degree 2, level 1 and weight k. By adelization (as described in more detail in the 
next section), there is a cuspidal automorphic representation irp canonically attached to F; the 
assumption that the level is 1 and there is no nebentypus means that ttf is unramified at finite 
places with trivial central character, as above. Thus we have local components tt p (F) £ X p and 
corresponding parameters (a p , b p ) G Y p for every prime p. 

The generalized Ramanujan conjecture has been proved in this setting by Weissauer |50j : it 
states that, if F is not a Saito-Kurokawa lift (these forms are defined in |13| for example; at the 
beginning of Section \5.2\ we recall the description in terms of L-functions) , we have tv p (F) G X for 
all p, i.e., \a p \ = \b p \ = 1. On the other hand, if F is a Saito-Kurokawa lift, then \a\ = 1 and 

{|6|,|6|- 1 } = {// 2 ,p- 1 / 2 }. 

Remark 1.4. Partly because our paper is meant to explore the general philosophy of families of cusp 
forms, we will not hesitate to use this very deep result of Weissauer when this helps in simplifying 
our arguments. But it will be seen that the proof of the local equidistribution property itself does 
not invoke this result, and it seems quite likely that, with some additional work, it could be avoided 
in most, if not all, of the applications (in similar questions of local equidistribution for classical 
Maass cusp forms on GL(2), one can avoid the unproved Ramanujan-Petersson conjecture). 

Denote by $* k any fixed Hecke-basis of Although this is not known to be unique, the aver- 
ages we are going to consider turn out to be independent of this choice. In fact, all final results 
could be phrased directly in terms of automorphic representations, avoiding such a choice (at least 
seemingly) . 

We next proceed to define our way of weighting the cusp forms in §>* k . This generalizes the 
ui k in the statement of the first theorem, and the reader may assume below that the parameters 
introduced are d = 4 and A = 1. 

Let d > be a positive integer such that — d is a fundamental discriminant of an imaginary 
quadratic field (i.e., one of the following holds: (1) d is congruent to 3 (mod 4) and is square-free; 
or (2) d = Am where m is congruent to 1 or 2 (mod 4) and m is square- free) . Let CLj denote the 
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ideal class group of this field, h(—d) denote the class number, i.e., the cardinality of CLj, and w(—d) 
denote the number of roots of unity. Finally, fix a character A of CLj. 

There is a well-known natural isomorphism, to be recalled more precisely in Section \2. 2\ between 
CLj and the SL(2, Z)-equivalence classes of primitive semi-integral two by two positive definite 
matrix with determinant equal to d/A. By abuse of notation, we will also use Chi to denote the set 
of equivalence classes of such matrices. 

Define normalizing factors 



c k,d 

where 



3 



dV 4c fc 



We note that, using Dirichlet's class number formula, one can also write 



^ = ^(4^-^-1)^-2). 



CkA 



cA 1 k 4-7TQ; 



w{-dfL{l, X d) 



where Xd is the real primitive Dirichlet character associated to the extension Q(V —d). Now, for 
each F G S£, we have a Fourier expansion 

F(Z) = ^a(F,T)e(Tr(TZ)), 

T>0 

where T runs over positive definite symmetric semi-integral matrices of size 2: 

a 6/2 N 



6/2 

with (a, b, c) G Z. It follows easily from the modularity property that a(F,T) depends only on 
the equivalence class of T modulo SL(2,Z). In fact, when k is even, a(F,T) depends only on the 
equivalence class of T modulo GL(2,Z). 
We now let 

F , \a(d,A;F)\ 2 

w fc,d,A - °k,d ■ "A • 



(F,F) 
where we put 

' ifA 2 = l, 



d A 
and 



otherwise. 



(1.2.3) a(d,A;F) = ^ A(c)a{F,c), 

cec\ d 

a quantity which is well-defined in view of the invariance of Fourier coefficients under SL(2,Z). 

Remark 1.5. We will often consider (d,A) to be fixed, and simplify the notation by writing cjjf 
only. Note that if d = 4 and A = 1, the weight ui^ is the same as the one defined in Theorem ll.il 

Next, we define the local spectral measures associated to the family S£; we will show that they 
become equidistributed as k — > +oo over even integers. Let S be a finite set of primes. We have 
the components tts(F) = (ir p (F)) pe s G X§, and we define the measure us,k on by 



dvs,k — dvs,k,d,A — ^2 ^k,d,A^ s {F) 



where 5, is the Dirac mass at the given point. As we will see, the normalization used has the effect 
that vs,k is asymptotically a probability measure on X$'- we will show that 

lim vs,2k(X s ) = 1. 

k— »+oo 

The local equidistribution problem — which can clearly be phrased for very general families of 
cusp forms — is to determine if these measures have limits as k — > +00, to identify their limit, to 
see in particular if the limit for a given S is the product of the limits for the subsets {p}, p G S 
(corresponding to independence of the restrictions), and finally — if possible — to express the 
resulting equidistribution in quantitative terms. 

To state our theorem, we now define the limiting measures. First of all, we define a generalized 
Sato- Tate measure /i on each X p by first taking the probability Haar measure on the space of con- 
jugacy classes of the compact unitary symplectic group USp(4), then pushing this to a probability 
measure on X by means of the map 



,i0i 



\ 



-i6i 



\ 



-i02 



^ (e id \e ie2 ) 



eX, (9\, 6 2 ) G [0, 7r] 2 , e 1 <9 2 , 



and finally extending it to X p by defining it equal to outside X. 

In terms of the coordinates (9i,9 2 ) on X, the resulting measure [x is explicitly given by 

(1.2.4) dn{e x ,Q 2 ) = -^(cos 6>i - cos 8 2 f sin 2 6>i sin 2 9 2 d9 1 d9 2 



IT" 



from the Weyl integration formula [13]. We can also interpret this measure as coming in the same 
way from conjugacy classes of the unitary spin group USpin(5,C), because of the "exceptional 
isomorphism" Sp(4) ~ Spin(5). 

For each finite set of primes S, and d, A as above, we now define the measure [is = ^s,d,A on 
Xs by the formula 

dfJ-s = dfi pAA 
pes 

where, for a single prime, we have 

dfj> p ,d,A = (1- (^)^a; 1 ^ 

(recall that /x is defined on X p , but has support on X only; the same is therefore true of /is) and 



the density functions A p = A 



(1.2.5) 



A,, 







((1 



'1 



P ,d,A are given by 

4 cos 2 1 ^ /^j _|_ 1 y, _ 4 cos 2 2 ' 

+ I(2 P cos0 lv ^-A p )( 2 ^ 



f((l + |) 2 - i£2 p)((l + |) 2 - i ^) 
1)2+ I(2ms0i./«- AJlf 2 ^! 

i) 2 + i(2cos0 2V p-A p )( 2 ^-A p )) 

2\ p cos 6>i — 2A p cos 02 , 1\ 



X 



((1 



if p inert, 

if p split, 
if p ramified, 



p/ \ vp P y 

where the behavior of primes refers to the field Q(\/—d), and in the second and third cases, we put 

A p= 2 A ( p )' 

7V(p)=p 

the sum over the (one or two) prime ideals of norm p in Q(\/— cZ). 

Although we have written down this concrete, but unenlightening, expression, there is a more 
intrinsic definition of these measures fj, p = fi Pt d,A and this will in fact be the way they will naturally 
occur (and the way we will prove the results): they are precisely what Furusawa and Shalika [18] call 
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the Plancherel measure for the local Bessel model associated to the data (d,A). In particular, this 
description shows that they are probability measures, which is not quite obvious from the definition. 
On the other hand, the following property, which is of great relevance to global applications, is 
immediate: as p — > +00, the measures fj, p converge weakly to the measure which has a group- 
theoretic interpretation. 

Our local equidistribution result can now be stated: 

Theorem 1.6 (Local equidistribution and independence). Fix any d,A as above. For any finite 
set of primes S, the measures v$,k on X$ converge weak-* to [is as k — ^ +00 over even integers, 
i.e., for any continuous function <p on Ys, we have 



lim V uj2 k (p{(a p (F),b p (F)) peS ) = / cp{x)dfi s (x) 

In particular, if 



pes 

is a product function, we have 



lim V u% k (p({a p {F),b p (F)) peS ) = TT / cp p (x)dfi p {x) 



FeS* pes 



Moreover, assume (p is of product form, and that (p p is a Laurent polynomial in (a, 6, a -1 , 
invariant under the action of the group W given by M.2.2}) . and of total degree d p as a polynomial 
in (a + a , b + b^ 1 ), then we have 



(1.2.6) V LU^((a p (F),b p (F)) peS )= I <p(x)dns(x) + 0(V 2 / 3 L 1+£ |M 

Fe$* J Ys v 

for any e > 0, where 



L=np d ?, 

pes 

and ||</?||oo is the maximum of \<p\ on X^ s \ C Ys- The implied constant depends only on d and e. 

Remark 1.7. It is possible to extend our results to odd k. However, this requires a slightly different 
definition of the weights u k . For simplicity, we only consider k even in this paper. 

Remark 1.8. In a recent preprint, Sug Woo Shin [37] has proved a related result. In Shin's work, 
the weights ui k are not present; instead the cusp forms are counted with the natural weight 1. 
Using the trace formula, he proves a qualitative result that for suitable families of cusp forms on 
connected reductive groups over totally real fields, there is local equidistribution at a given place; 
when the level grows the sum of the point measures associated to the forms of fixed level converges 
towards the Plancherel measure on the unitary dual of the local group. 

In fact, from the viewpoint of automorphic representations on reductive groups, our result is 
essentially a (quantitative) relative trace formula analogue of what Shin (and others before him, 
such as DeGeorge-Wallach [9], Clozel jSJ, Savin [35], Serre |46j . Sauvageot [H]) did using the trace 
formula. 

We expect that the methods of this paper would suffice to prove the equidistribution result for 
Siegel modular forms of level N coprime to the set of places S as N + k — > 00. It would also 
be interesting to see if our results can be generalized to the case of automorphic forms on the 
split special orthogonal groups, using the formulas for the Bessel model there from [6] (at least 
qualitatively). We hope to treat these questions elsewhere. 



We briefly explain the structure of this paper. In Chapter [21 we introduce the Bessel model, 
explain its relation to the Fourier coefficients and derive a result relating the Fourier coefficients to 
Satake parameters. In Chapter[3]we recall the definition of Poincare series in this context and derive 
a Petersson-type quantitative orthogonality formula for the Siegel cusp forms (this involves non- 
trivial adaptations of the method of Kitaoka [26J). In Chapter HI we put the above results together 
to deduce our main theorem (Theorem II. 6p on local equidistribution. Finally, in Chapter [5l we 
prove Theorems ll.ll and ll.2l as well as provide several other applications of the results of the previous 
chapters. 

Acknowledgements. We would like to thank Masaaki Furusawa for forwarding us a copy of 
the relevant part of his ongoing work [T7j, and for some helpful suggestions. Thanks also to a 
referee for prompting us to give a more thorough explanation of the behavior of Theorem 11.21 

1.3. Notation. We introduce here some notation used in the paper. 

• The symbols Z, Z>o, Q, R, C, Z p and Q p have the usual meanings. A denotes the ring of 
adeles of Q. For a complex number z, e(z) denotes e 2mz . 

• For any commutative ring R and positive integer n, M(n, R) denotes the ring of n by n 
matrices with entries in R and GL(n, R) denotes the group of invertible matrices in M(n, R). 
If A G M(n, R), we let 1 A denote its transpose. We use R x to denote GL(1, R). 

• For matrices A and B, we use A[B] to denote l BAB, whenever the matrices are of compatible 
sizes. 

• We say that a symmetric matrix in M(n, Z) is semi-integral if it has integral diagonal entries 
and half-integral off-diagonal ones. 

• Denote by J n the 2n by 2n matrix given by 

o V 



We use J to denote J2. 

• For a positive integer n, define the algebraic group GSp(2n) over Z by 

GSp(2n,R) = {g£ GL(2n, R)\ l gJ n g = ii n {g)J n ^ n {g) G R x } 

for any commutative ring R. 

Define Sp(2n) to be the subgroup of GSp(2n) consisting of elements g\ G GSp(2n) with 

t*n(gi) = i- 

The letter G will always stand for GSp(4). The letter T will always stand for the group 
Sp(4,Z). 

• The Siegel upper-half space is defined by 

M n = {Z G M n (C)\Z = l Z, Im(Z) is positive definite}. 



For 



and Z G Efe, we denote 



J{g,Z) = CZ + D. 

• For a prime p, the maximal compact subgroup K p of G(Q P ) is defined by 

K p = G(Q p )nGL(4,Z p ). 

• For a quadratic extension % of Q and p a prime, define % p = % 0q Q p ; we let Z^ denote 
the ring of integers of % and ^x,p its p-closure in % p . 

9 



2. Bessel models 



2.1. Global Bessel models. We recall the definition of the Bessel model of Novodvorsky and 
Piatetski-Shapiro |35j following the exposition of Furusawa |15j . 

Let S G M2(Q) be a symmetric matrixQ Define disc(S) = — 4det(S) and put d = — disc(5). For 



(a b/2\ 
S -{b/2 c) 



we define the element 

5 = &-'' V2 c 



-a -b/2 y 

Let % denote the subfield Q(V~ d) of C. We always identify Q(£) with % via 

(2.1.1) Q(0 3x + y^x + y^^ G 3C, x,y G Q. 

We define a Q-algebraic subgroup T = T$ of GL(2) by 

(2.1.2) T = {g G GL(2)|^ = det(g)S}, 

so that it is not hard to verify that T(Q) = Q(£) x - We identify T(Q) with 3C X via (j2"XTj) . 
We can also consider T as a subgroup of G via 

Let us further denote by U the subgroup of G defined by 

U = {u(X) = ^ f)\ t X = X}, 

and by R be subgroup R = TU of G. 

Let ^ = ip v be a character of A such that 

• The conductor of ip p is Z p for all (finite) primes p, 

• i>oo(x) = e(x), for 

• = 1. 

We define the character 6 = 6s on U (A) by 

0(u(X)) = i/,(Tr(S(X))). 

Let A be a character of T(A)/T(Q) such that A|A X = 1. Via (|2.1.ip we can think of A as a 
character of D<C X (A)/3<; X such that A|A X = 1. Denote by A (g> 9 the character of i?(A) defined by 
(A (8) = A(t)0(it) for i G T(A) and u G t/"(A). 

Let 7T be an automorphic cuspidal representation of G(A) with trivial central character and V w 
be its space of automorphic forms. For $ G 14- , we define a function 5$ on G(A) by 



(2.1.4) B*(h)= (A®0)(r)*(rfc)dr. 



The C- vector space of functions on G(A) spanned by {-B$| ( 3? G V^-} is called the global Bessel 
space of type (S, A, ip) for ir; it is invariant under the regular action of G(A) and when the space 
is non-zero, the corresponding representation is a model of 7r. Thus one says that tt has a global 
Bessel model of type (S, A, ^) if this global Bessel space is non-zero, i.e., if there exists $ G V n such 
that B§ ^ 0. 



^ The notation conflicts a bit with the sets of primes S, but we hope that the boldface font of the latter and the 
context will prevent any confusion. 
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r*(7V) :={g € T:g = 



(mod TV)}. 



2.2. The classical interpretation. Let us now suppose that tt come from a classical Siegel 
cusp form F. More precisely, for a positive integer N, define 

/* o\ 
0*00 
0*0 
\0 */ 

We say that F G Sfc(r*(iV)) if it is a holomorphic function on H2 which satisfies 

F( 7 Z) = det(J( 7 ,Z)) fc F(Z) 
for 7 G r*(iV), Z G H2, and vanishes at the cusps. It is well-known that F has a Fourier expansion 

F(Z)= a(F,T)e(Tv(TZ)), 

T>0,Te£ 

where e(z) = exp(2Triz) and T runs through all symmetric positive-definite matrices of size two in 
a suitable lattice L (depending on N). If TV = 1, then L is just the set of symmetric, semi-integral 
matrices. Also, recall that Sfc(r*(l)) is denoted simply by 

We define the adelization $p of F to be the function on G(A) defined by 

(2.2.1) 3> F {lhM = / u 2 (/ ioo ) fc det(J(/ ioo ,i/ 2 ))- A; F(/ loo (i)) 

where 7 G G(Q),/ioc G G(R)+ and 

p\N p\N 

where 

/* 0\ 



< = {g G G(Zp) : 5 = 



(mod N)} 



0*00 
0*0 
\0 */ 

is the local analogue of F*(N). Then is a well-defined function on the whole of G(A) by strong 
approximation, and is an automorphic form. 

We now assume TV = 1. Let d be a positive integer such that —d is the discriminant of the 
imaginary quadratic field Q(\/— ~d). Define 



(2.2.2) 



S = S{-d) = < 



t 0" 
1, 



l+rf 1 
4 2 

\ \ 



if d = (mod 4), 
if d ee 3 (mod 4). 



Define the groups i?, T, J7 as in the previous section. 

Put % = Q(y/—d). Recall that CLj denotes the ideal class group of X. Let (t c ), c G CLj, be coset 
representatives such that 



(2.2.3) 

with t c G n p <oo 

with 7c G GL(2,( 
The matrices 



T(A) = JJt c T(Q)T(IR)n p<00 (T(Q p ) n GL 2 (Zp)), 

c 

. We can write 

t c = ^j c m c K c 
m c G GL + (2,M), and k c G n p<00 GL(2, Z p ). 

S c = det( 7c )" 1 ('7c)57c 
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have discriminant —d, and form a set of representatives of the SL(2, Z)-equivalence classes of prim- 
itive semi-integral positive definite matrices of discriminant —d. 

Choose A a character of T(A)/T(Q)T(R)((U P<00 T(Z P )), which we identify with an ideal class 
character of <Q(V—d). 

Next, for any positive integer N, define (this is a certain ray class group) 



C\ d (N) = T(A)/T(Q)T(M)n p<00 (T(Q p ) n K^(N)), 



where K p °\n) is the subgroup of GL2(Z p ) consisting of elements that are congruent to a matrix 
of the form 

'* 0\ , , 

(mod N). 



* * 

As before, we can take coset representatives (i c >), d G CLj(AT), such that 

(2.2.4) T(A) = Y[t c/ T(q)T(R)U p<0O (T(q p )nK^(N)), 

d 

with t d G n p <oo T ( ( Qp)' and write 

t c , = -y c/ m c >K c / 

with 7c / G GL(2,Q), m c / G GL+(2,M), and k c > G n p<oo ^ 0) (A^). 
Define the matrices 

S c > = det( 7c _1 (*7c')^7c- 
For each pair of positive integers L, M, define the element H(L, M) G G{A) by 

(LM 2 \ 
LM 



H(L, M)oo = 1, H(L, M) p 



1 



for each prime p < oo. Note that H(l, 1) = 1. 

For any symmetric matrix T, we let T L,M denote the matrix 



(2.2.5) T^ = ^ JT^ y 
Define the quantity B(L, M;® F ) by 

(2.2.6) B(L,M;$ F ) = B*(H(L,M)) = [ (A <g> 0)(r)$ F (rH(L, M))dr. 

JR(Q)Z G (A)\R(A) 

The next proposition proves an important relation which expresses the quantity B(L, M;<& F ) 
in terms of Fourier coefficients. The proof is fairly routine, but to our best knowledge it has not 
appeared in print before. 

Proposition 2.1. Let F G have the Fourier expansion 

F(Z) = J2a(F,T)e(Tr(TZ)). 

T>0 

Then we have 

B(L,M^ F ) = r.e~^ s \LM)- k —^-— £ A^)a(F,^ M ). 

where r is a non-zero constant depending only on the normalization used for the Haar measure on 
R. 
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Proof. For the purpose of this proof, we shorten H(L, M) to H whenever convenient. Note that 

B(L,M;<f> F ) =B(1,1;^' M ) 
where the automorphic form <& F ' M is given by 

$ L /\g) = $ F (gH). 



Define 

(LM 2 \ 
LM 

1 



Hoc 



g G(M.y 



and define the Siegel modular form 

(2.2.7) F'{Z) = (LM)- k F(H- l Z), 

which is in S/c(r* (LM 2 )), as one can easily show. 

Let $f' be the adelization of F' as defined by (|2,2.ip . We claim that & F ' = &f'- To see this, 
since both functions are right invariant under the group 

n k p n 



LM 2 
V 

p\LM 2 p\LM 2 



it is enough to show that ^ F ,M (g OQ ) = &f'(9oo) for #oo G This is shown by the following 

computation: 

® L F M (9oo) = $ F (9ooH) 

= &F(Hn9°o) 

= f i 2 (H^ 1 g O0 ) k det(J(H^ 1 g O0 ,iI 2 )r k F(H^g O0 (i)) 
= (LM)~ k ^{g 00 ) k det{J{g 00 ,ih))- k F(H^ l g 00 (i)) 

= $F'(5oc). 

Hence we are left with the problem of evaluating 5(1, 1; &f')- Note that <&p' is right invariant 
under K^\m) (where we think of GL2 as a subgroup of GSp 4 via (|2.1,3p ). Using (|2.2.4p . the same 
arguments as in \42\ Prop. 2.8.5], give us 

B(l,l;^) = e- a ^( g ) 1 Yl W)a(F',S c ) 

1 d[ n ceCi d (M) 

for a suitably normalized Haar measure, where a(F',T) denotes the Fourier coefficients of F' . 
Using (|2.2.7p . one can easily check that 

a(F',S c ) = (LM)- k a(F,St M ), 

and this completes the proof. □ 

Remark 2.2. In an earlier preprint version of this paper, we had claimed such a result with a sum 
over c in CLj instead of CLj(M). This was incorrect (when M 7^ 1); the mistake in the proof was 
to assume that $f' 1S invariant under the bigger subgroup GL2(Z p ) when arguing as in \42\ Prop. 
2.8.5]. 

Remark 2.3. The above result is one of the three crucial ingredients that are required for the 
proof of the asymptotic Petersson-type formula (Proposition 13. 6p which forms the technical heart 
of this paper. The other two ingredients are Sugano's formula (Theorem I2.5P and the asymptotic 
orthogonality for Poincare series (Proposition 13. 3| . 
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2.3. Local Bessel models and Sugano's formula. Let it = ® v tt v be an irreducible automorphic 
cuspidal representation of G(A) with trivial central character and V v be its space of automorphic 
forms. We assume that tt is unramified at all finite places. Let S be a positive definite, symmetric, 
semi-integral matrix such that — d = —4det(S) is the discriminant of the imaginary quadratic field 
DC = Q(y/—d). Let ip, A be defined as in the previous Section. Define the groups R, T, U as before 
and the Bessel function B§ on G(A) as in (|2.1.4p . for a function $ = Y\ v <& v which is a pure tensor 
in 7r. 



For a finite prime p, we use y—J to denote the Legendre symbol; thus y—J equals —1, or 1 
depending on whether the prime is inert, ramified or split in %. In the latter two cases, we use p% 
to denote an element of % p = % (8>q Q p such that Nx/q(px) £ P^p ■ 

Outside v = oo, the local representations are unramified spherical principal series. Therefore, by 
the uniqueness of the Bessel model for G, due to Novodvorsky and Piatetski-Shapiro [35J, we have 

(2.3.1) B*(g) = B^goc) J] B p(9p) 

p 

where B p is a local Bessel function on G(Q P ), the definition of which we will now recall. 

Remark 2.4. If the global Bessel space is zero, then both sides of (I2.3.ip are zero. In particu- 
lar, (|2.3.ip remains valid regardless of whether our choice of S and A ensures a non-zero Bessel 
model. 

To describe the local Bessel function B p for a prime p, let S be the space of locally constant 
functions (p on G(Q P ) satisfying 

<p(tuh) = A p (t)9 p (u)ip(h), for t G T(Q p ),u G U(Q p ),h G G(Q P ). 

Then by Novodvorsky and Piatetski-Shapiro [35], there exists a unique subspace 23(7r p ) of ¥> such 
that the right regular representation of G(Q P ) on 23(7r p ) is isomorphic to tt p . Let B p be the unique 
-fTp-fixed vector in 23(-7r„) such that B p (l) = 1. Therefore we have 

(2.3.2) B p (tuhk) = A p {t)6p(u)B p (h), 

for t G T(Q p ),u G U(Q p ),h G G(Q P ), fc G i^. 

Let h p (l,m) G G(Q P ) be the matrix defined as follows: 

h p (l, m) : = 



pl+m 



1 

\ p m J 

As explained in |15| . the local Bessel function B p is completely determined by its values on 
h p (l,m). An explicit formula for B p (h p (l,m)) in terms of the Satake parameters is stated in [6]. 
This formula can be neatly encapsulated in a generating function, due to Sugano |48j . which we 
now explain. 

Because ir p is spherical, as recalled earlier, it is the unramified constituent of a representation 
Xi x X2 * o induced from a character of the Borel subgroup associated to unramified characters 
Xi) X2, v of Q p ) and because it has trivial central character (since tt does) we have X1X2C 2 = 1- Let 
us put (a p , bp) = (a(p),a(p)xi(p)) G Y p , and as in the definition of the measure fi p , let 



where the number of terms in the sum is 1 + , 



N(x)=p 
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The next Theorem is due to Sugano [Ml p. 544] (the reader may also consult [15j (3-6)]). 

Theorem 2.5 (Sugano). Let tt be an unramified spherical principal series representation ofG(Q p ) 
with associated local parameters (a, b) € Y p and spherical Bessel function B p as above. Then we 
have 

(2.3.3) B p (h p (l,m))=p~ 2m -%U l p ' m (a,b) 
where for each I, m ^ 0, the function 

U^ m (a,b) = U^ m (a,b;X p ,A p ) 
is a Laurent polynomial in C[a,b, a ,b ], invariant under the action of the Weyl group hi. 2.2(1 , 
which depends only on p, & n d \>- 

More precisely, the generating function 

(2.3.4) C P (X, Y) = C P (X, Y;a,b) = ^Yl U p' m ( a > 

l>0 TO>0 

is a rational function given by 
where 

P P (X) = (1 - abX)(l - ab~ l X){l - a _1 6X)(l - a^b^X), 
Q P (Y) = (1 - Y)(1 - bY){\ - a- 1 y)(l - b~ x Y), 
H P (X,Y) = (1 + XY 2 )(M 1 (X)(1 + X) +p- 1 / 2 X p a(a,b)X 2 ) 

-XY(a(a,b)Mx(X) - p- 1/2 X p M 2 (X)) - p^ 2 X p P p (X)Y + p' 1 P P (X)Y 2 , 

in terms of auxiliary polynomials given by 

a(a, b) = a + b + a" 1 + b' 1 , r(a, b) = 1 + ab + ab' 1 + oT l b + a~ x b~ x , 



M X {X) = 1- [p 



1K X x ~ ] 



P V 2 \ p o-(a, b) - (r(o, b) - 1) - \£j X - p- 1 (^j X 2 , 



P . 

M 2 (X) = 1 - r(a, b)X - r(a, b)X 2 + X 3 . 
Remark 2.6. For instance, we note that 

£/°'°M) = l 

and that 

(2.3.6) U x '°{a, b) = a (a, b) - p~ 1/2 \ p = a + b + a' 1 + - p~ 1/2 \ p . 

We also note that taking X = leads to the simple formula 



i 



l-p-i\ p Y +p-i(|)r 2 



(2.3.7) E^(^)V Qp(y) 

which we will use later on. The formula for Y = is more complicated, but we note (also for 
further reference) that it implies the formula 

(2.3.8) U^(a,b) = r(o,6) - (p - ' (p X ' 2 X p a{a,b) - (j) (r(a,b) - 1) - . 
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As was the case for the definition of the measures fx p , we have written down a concrete formula. 
These are not very enlightening by themselves (though we will use the special cases above), and 
the intrinsic point of view is that of Macdonald polynomials [32] associated to a root system. In 
particular, this leads to the following important fact: 

Proposition 2.7. Let (d,A) be as before. For any fixed prime p, the functions 

(a,b)^U l p m (a,b) 

where I, m run over non-negative integers, form a basis of the space of Laurent polynomials in 
C[a, b, a" 1 , 6 _1 ] which are invariant under the group W generated by the three transformations 
above. 

Moreover, any such Laurent polynomial ip which has total degree d as polynomial in the variables 
(a + a -1 , b + b^ 1 ) can be represented as a combination of polynomials Up" 1 (a, b) with I + 2m ^ d. 

Proof. Because of Theorem 12.51 we can work with the Bessel functions B p (h p (l,m)) instead. But 
as shown in detail in [17, §3], these unramified Bessel functions are (specializations of) Macdonald 
polynomials associated to the root system of G, in the sense of [32] . By the theory of Macdonald 
polynomials, these unramified Bessel functions form a basis for the Laurent polynomials in two 
variables that are symmetric under the action of the Weyl group W. 

The last statement, concerning the U p m occurring in the decomposition of ip of bidegree (d,d), 
can be easily proved by induction from the corresponding fact for the coefficients (say U^ m {a,b)) 
of the simpler generating series 



for which the stated property is quite clear. (It is also a standard fact about the characters of 
representations of USp(4, C), since a(a,b) and r(a,b) are the characters of the two fundamental 
representations acting on a maximal torus.) □ 

Lemma 2.8. Let (d,A) be as before. Let S be a finite set of primes and (l p ), (m p ) be tuples of 
non-negative integers, indexed by S. There exists an absolute constant C ^ such that for every 
(xp)pgs = (a p ,b p ) G X s , i.e., parameters of tempered representations, we have 

| H U^iap, ftp) | < J] (l p + 3) 3 (m p + 3) 3 . 
pes pes 



Proof. It is enough to prove this when S = {p} is a single prime, and l p = I, m p = m ^ 0. Then 
by Sugano's formula, the polynomial U p m (a,b) is a linear combination of at most 14 polynomials 
of the type arising in the expansion of the denominator only, i.e., of 



P p (X)Qp{YY 

and moreover the coefficients in this combination are absolutely bounded as p varies (they are either 
constants or involve quantities like p~ l l 2 ). 

Expanding in geometric series and using \a\ = 1, \b\ = 1, the coefficient of X m Y l in the expansion 
of the denominator is a product of the coefficient of X m and that of Y l ; each of them is a sum, with 
coefficient +1, of ^ (m + 3) 3 (resp. (/ + 3) 3 ) terms of size ^ 1. The result follows from this. □ 

Remark 2.9. Sugano's formula explicitly computes the Bessel function in terms of Satake parameters 
in the case of an unramified representation. The other case where an explicit formula for the Bessel 
function at a finite place is known is when tt p is Steinberg, see [41], [37] . 
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2.4. The key relation. We consider now Siegel modular forms again. Let 

F(Z) = a(F,T)e(Tr(TZ)) G § k 

T>0 

be an eigenfunction for all the Hecke operators. Define its adelization <&f(9) by (|2.2.ip . This is a 
function on G(Q)\G(A) and we may consider the representation of G(A) generated by it under the 
right-regular action. Because we do not have strong multiplicity one for G, we can only say that 
this representation is a multiple of an irreducible representation irp. However, the unicity of ttf, 
as an isomorphism class of representations of G(A), is enough for our purposes!! 
We can factor irp = <gnr v (F) where the local representations ir v are given by: 



holomorphic discrete series if v = oo, 

unramified spherical principal series if v is finite, 



and we denote by (a p (F) , b p (F)) € Y p , the local parameters corresponding to the local representa- 
tion TTp(F) at a finite place. 

Let once more d be a positive integer such that —d is a fundamental discriminant and define S 
as in (|2.2.2p . Choose an ideal class character A of DC. Let the additive character ip, the groups 
R, T, U and the matrices S c , S c ,' be defined as in Section 12.21 For positive integers L, M, define 
B(L,M; &f) by (|2.2.6I) . Then, by the uniqueness of the Bessel model (i.e., (|2.3.ip . we have 

B(L, M; $>f) = B(l, 1; $f) JJ B p (h p (l p , m p )) = B(l, 1; <f> F ) H B p (h p (l p ,m p )), 

V p\LM 

where L and m p are the p-adic valuations of L and M respectively. Now, using Sugano's for- 
mula (]2.3.3p and twice Proposition 12.11 - which has the effect of canceling the constant r ^ that 
appears in the latter -, we deduce: 

Theorem 2.10. Let (d,A) be as before, let p be prime and let U p ' m (a,b) be the functions defined 
in Theorem \2.5[ For any F E §2 fc and integers L, M ^ 1, we have 



A(d)a(F,S^ M ) = L k -lM k - 2 £ A(c)a(F,S c ) ]J U 1 /'^ (a p (F),b p (F)), 



I CUM) I ^ 

1 K n c'eCl d (M) ceCl d p\LM 

where l p and m p are the p-adic valuations of L and M respectively. 

The point of this key result is that it allows us to study functions of the Satake parameters of 
Tip using Fourier coefficients of F, although there is no direct identification of Hecke eigenvalues 
with Fourier coefficients. 

Remark 2.11. This relation holds for every A, but we can not remove the sum over c by Fourier 
inversion because the functions U p " mp depend on A. 

3. POINCARE SERIES, PETERSSON FORMULA AND ORTHOGONALITY 

The relation given by Theorem 12.101 between Fourier coefficients of F on the one hand, and func- 
tions of the spectral Satake parameters of irp on the other, will enable us to deduce equidistribution 
results for Satake parameters from asymptotics for Fourier coefficients. For this, we need a way 
to understand averages of Fourier coefficients of Siegel forms F in a suitable family; this will be 
provided by a variant of the classical Petersson formula. In order to prove the latter, we follow the 
standard approach: we consider Poincare series and study their Fourier coefficients. 



Added in proof: in a recent preprint, Narita, Pitale and Schmidt show that $f does indeed generate an irreducible 
representation. 
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3.1. Poincare series and the Petersson formula. Given a symmetric semi-integral positive- 
definite matrix Q of size two, the Q-th Poincare series of weight k, denoted P k) Q, is defined as 
follows: 

P k , Q {Z)= Y, det(J( 7 ,Z))- fc e(Tr(g 7 (Z))) 

7£A\r 

where A is the subgroup of T consisting of matrices of the form ( \ 1 ) , with U symmetric. 



l / 

It is known that P k ,Q is absolutely and locally uniformly convergent for k 6, and defines 
an element of §> k (as first proved by Maass). In fact, any Siegel cusp form F G § k is a linear 
combination of various Pk,Q (with Q varying). This follows from the basic property of Poincare 
series: they represent, in terms of the Petersson inner product, the linear forms on § k given by 
Fourier coefficients. Precisely, for F E §>k with Fourier expansion 



F(Z) = J>(F,T)e(Tr(TZ)), 



T>0 

we have the crucial identity 

(3.1.1) (F,Pk,T) = 8c k (detT)- k+3 / 2 a(F,T), 
where 

(3.1.2) c k = ^{^f- 2k T{k - l)Y{k - 2). 

(see [28J or [27, p. 90] for instance). 

We are interested in the limiting behavior of a(k; c,d ,L, M) as k — > +oo. The following qualita- 
tive result was proved in [31] : 

Proposition 3.1 (Asymptotic orthogonality, qualitative version). For L, M ) 1, c £ Cl^ and 

d £ Cl d (M), let 

a(k;c,c ,L,M) = a(P k , Sc , S^ ,M ) 
denote the S^,' M -th Fourier coefficient of the Poincare series P k ,s c - Then, we have 

a(k;c,d,L,M) -> | Aut(c)| • 5(c,c;L;M) 
as k — )• +oo over the even integers. Here 

1 if L = 1,M = 1 and c is GL(2,Z) -equivalent to d, 



5(c,c';L;M) 



otherwise, 



and | Aut(c)| is the finite group of integral points in the orthogonal group 0(T) of the quadratic 
form defined by c. 

Remark 3.2. Since this will be a subtle point later on, we emphasize that 5(c, d; L, M) = 1 when c 
and d are invariant under GL(2, Z), not under SL(2,Z). 

This is sufficient for some basic applications, but (for example) to handle the low-lying zeros, we 
require a quantitative version. We prove the following: 

Proposition 3.3 (Asymptotic orthogonality, quantitative version). With notation as above, we 
have 

a{k;c,c',L,M) = |Aut(e)| • 5{c,d;L;M) + L k - 3/2 M k - 2 A(k;c,c ,L,M) 

where 

A{k;c,c',L,M) < L 1+e M 3 / 2+e £T 2 / 3 
for any e > 0, the implied constant depending only on £ and d. 
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In the proof, for conciseness, we will write \A\ for the determinant of a matrix. The basic 
framework of the argument is contained in the work of Kitaoka [26], who proved an estimate 
for the Fourier coefficients a(Pk t Q,T) of Poincare series for fixed Q and k 6, in terms of the 
determinant det(T) (and deduced from this an estimate for Fourier coefficients of arbitrary Siegel 
cusp forms in since the space is spanned by Poincare series). 

However, Kitaoka considered k to be fixed; our goal is to have a uniform estimate in terms of L, 
M and k, and this requires more detailed arguments. 

In particular, we will require the following quite standard asymptotics for Bessel functions: 



(3.1.3) J k {x) < — if k > 1, 0«Cx« Vk + 1, 

T{k + 1) 

(3.1.4) J k (x) <min(l,x/r 1 )Ar 1 / 3 , if k ^ 1, x > 1, 

(3.1.5) Jfc(x) < iffc^l, x>0, 



where the implied constants are absolute (the first inequality follows from the Taylor expansion 
of Jk(z) at z = 0, the second is [24J, (2.11)], and the third, which is very rough, by combining 
\J k (x)\ < 1 when x < 2k, and, e.g., [Ml (2.11')] for x ^ 2fc). 

Proof of Proposition \3.3X Let 

T = S^,' M , Q = S c , 

so that we must consider the T-th Fourier coefficients of Pk,o (this notation, which clashes a bit 
with the earlier one for the torus T, is chosen to be the same as that in [26], in order to facilitate 
references). Before starting, we recall that since we consider d to be fixed, so is the number of ideal 
classes, and hence Q varies in a fixed finite set, and may therefore be considered to be fixed. Also 
note that 

(3.1.6) det(T) = dL 2 M 2 /4, 

and we seek estimates involving det(T). Thus, compared with Kitaoka, the main difference is the 
dependency on k, which we must keep track of. In particular, we modify and sharpen Kitaoka's 
method, so that any implicit constants that appear depend only on d. 

Because we think of d as fixed, throughout the proof we drop the subscript d from the symbols 
<C, S>, x. The reader should not be misled into thinking that the implied constants are independent 
of d. 

Since the proof is rather technical, the reader is encouraged to assume first that d = 4, M = 1 
(so that there is a single class c = c', and moreover Q = S c = S c i = 1) and alscll by A2.2.5 j> . T is a 
simple diagonal matrix 

' L N 
,° L s 

In principle, we now follow the formula for a(Pk,Q, T) which is implicit in [26]. Given a system of 
representatives f) of ri(oo)\r/ri(oo), Kitaoka defines certain incomplete Poincare series Hk(M,Z) 
such that 

Pk, Q (Z)= Y,H k (M,Z). 

M6I) 

Denoting the T-th Fourier coefficient of (M, Z) by hk (M, T) , we have 

a(k; c, c',L, M) = a{P kiQ ,T) = ^ h k (M, T). 

Mel) 



This is the only case needed in Theorem 11.11 for averaging the spin L-function; however, this is not sufficient for 
Theorem 11.21 although the latter is also concerned only with the spinor L-function. 
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We write 



M 



A B 
C D 



where A, B, C and D are matrices in M(2, Z), and we now divide the sum above depending on the 
rank of C. We denote the component corresponding to rank i by Ri, so that 

a(P k ,Q,T) =R + R 1+ R 2 . 

Step 1 (rank 0). First of all, we consider Rq. By [26\ p. 160] (or direct check), we have 



Ro 



£/eGL(2,Z) 
UT t U=Q 



which is unless T is GL(2, Z)-equivalent to Q, in which case it is equal to | Aut(T)| = \0(T, Z)| 
(where T is viewed as defining a quadratic form and 0(T) is the corresponding orthogonal group). 
In our case, looking at the determinant we find that Rq = unless L = M = 1, and then it is also 
except if c is GL(2, Z)-equivalent to c', and is then | Aut(c)|. In other words, we have 

R = | Aut(c)| ■ 6{c,c';L;M), 

and hence, by definition, the remainder is therefore 

(3.1.7) R l+ R 2 = L k -^ 2 M h - 2 A(k- c, c', L, M), 

and — having isolated our main term — we must now estimate the two remaining ones. 

Step 2 (rank 1). Following the computations in Kitaoka (specifically, Lemma 4, p. 159, Lemma 
1, p. 160, and up to line 2 on p. 163 in [2B]), but keeping track of the dependency on k by keeping 
the factor Q 3//4_fc / 2 (which Kitaoka considers as part of his implied constant), we find that 



(3.1. 



I#i|« £ E m^-^lQl^^Aim^m-^im^) 1 



/2 



■4-3/2 I 47T 



VWWh 

mc J 



where A(m, T) is the number of times T, seen as a quadratic form, represents m. 

Now recall that \Q\ = d/A and \T\ = L 2 M 2 (ci/4), and observe that A(m, T) = unless L divides 
m and A(m,T) = A(m/L, S c ) whenever L divides m. It follows that 

A(m,T) < £ {m/Lf 

for any e > 0. Using (|3.1.8p . we get by a very rough estimate that 



\Ri\ < £ (LM) k -l rn-^^im/Lfim^c) 1 / 2 



L\m 



J, 



k-3/2 



LMd\ 



Ti- 



me 



< {LM) k ~lL £ 



711 



-1/2+e 



mi, c 



,V2 



Jh 



Md\ 



fc— 3/2 ( 7T I 



Now we define 



m-^+Hm^c) 1 ' 2 



c,m^l 
= + ^12 + ^13 



J/ 



fc-3/2 



Md^ 



Ti- 



me / 
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where Jin corresponds to the sums restricted to 

mc > irMd, if i = 1 

■KMdk- 1 / 2 ^mc^ irMD, Hi = 2 
mc^TxMdk- 1 / 2 if i = 3. 

For i = 1 , the argument of the Bessel function is ^ 1 and by (|3.1.3[) , we find 



mOnM d 



We can replace the exponent k — 3/2 in the sum with any exponent 1 + 5, for small < 5 ^ 1 
(since k ^ 6 anyway), and then we can remove the summation condition, observing that the double 
series is then convergent, and obtain 

•Jin < M 1+£ k~ E , 

(taking 5 small enough in terms of e) for any E ^ 1 and e > 0, where the implied constant depends 
on E, e and d. 

For i = 2, we use (|3.1,3p and find that 

jLfc/2-3/4 

fti2 « * _ Y, rn-^(m,c)^«M^k- E 

^ ' ' cm^irMd 

for E ^ 1 and e > again (by summing over c first and then over to, and by Stirling's formula). 
Finally, using now (|3.1,4p . we have 

*13 «, k' 1 / 3 Y, ™- 1/2+£ (™, c) 1 ' 2 « k- l ^- l l 2 M l+ ' = k~ 5 / 6 M l+£ , 

cmK-irMdk- 1 / 2 

(summing as for ^12) • 
It follows that 

#i = R u + R 12 + 3?i 3 < M 1+£ A;-f 
for any e > 0, and so the contribution of rank 1 is bounded by 

(3.1.9) 122x1 <C (LM) fc ~§L E M 1+E Art 

for any e > 0, where the implied constant depends only on d and e. 

Step 3 (rank 2). Finally, we deal with the R 2 term, which is much more involved. The relevant 
set of matrices M is given by 



M G 



" *^}cSp(4, 



C D; 

where \C\ 7^ and D is arbitrary modulo C . Denoting 

M*(Z) = {CeM 2 (Z) I \C\?0}, 

we have then 

#2= E hk(M, T). 

CeMJ(Z) D (modC) 

The inner sum was computed by Kitaoka [26, p. 165, 166]. To state the formula, let 

p = p(C) := TQ^C- 1 } = T( l C- l )QC- 1 , 

and let 

(3.1.10) < sx s 2 
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be such that s 2 , s 2 are the eigenvalues of the positive definite matrix P. Then Kitaoka proved that 

(3.1.11) £ h k (M,T) = ^(g|) fc/2 " 3/ Vr 3/2 ir(Q,r;C)a fc (P(C)), 

where K(Q,T;C) is a type of matrix-argument Kloosterman sum (see [26, §1, p. 150] for the 
precise definition, which we do not need here), ancQ 

,7T/2 

3k(P) = / J A: _3/2(47rsisin6')J fe _3/2(47rs2sin6')sin6'd6'. 
Jo 

We note that 

\P\ = |T||Q||C|- 2 = (d/4) 2 L 2 M 2 \C\~ 2 . 

In order to exploit this formula (|3.1.1ip . we must handle the sum over C. For this purpose, 
we use a parametrization of M| (Z) in terms of principal divisors: any C G M| (Z) can be written 
uniquely 

(3.1.12) C=U ~ 1 ( C Q c 2 ) F_1 
where 

Uq, ci | ca, U £ GL(2,Z) and V £ SL(2,Z)/T° (<&/<*), 
where T°(n) denotes the congruence subgroup of SL(2, Z) (conjugate to To(n)) consisting of matrices 

, c d, 

with n | 6. Note that there is a bijection 

SL(2,Z)/r°(n) ~ P 1 (Z/nZ), 
(this is denoted 5(n) in [26]) and in particular 

(3.1.13) |SL(2,Z)/r°(n)| =n[](l+p- 1 ) <n 1+e 

p\n 

for any e > 0. 

We will first consider matrices where the last three parameters c = (c\,C2,V) are fixed, subject 
to the conditions above. The set of such triples is denoted V, and for each c £ V, we fix (as we 
can) a matrix U\ G GL(2, Z) such that the matrix 

A(c) = A:=T[v{^ 1 \ XJ X 
is Minkowski-reduced. This matrix is conjugate to a diagonal matrix H = H(c) of the form 

H=( a 

with a < c. Computing determinants and using the fact that A is Minkowski-reduced, we note also 
that we have 

(3.1.14) (d/4)__ =ocX8fai = (d/4) 2 — j-g- 

c l c 2 c l c 2 



(we recall again that d is assumed to be fixed). 



^We have made the change of variable t = sin(#) for convenience. 
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For fixed c S V, the set of matrices C E M| (Z) corresponding to c can be parameterized in the 
form 

where U varies freely over GL(2,Z) (this is a simple change of variable of the last parameter 
U £ GL(2,Z) in (13.1.12ft ). As shown in \26\ p. 167], for any such C associated to c, we have also 

|P| x \A\, Tr(P) x Tr(A[U}) = Tr(H[V]). 

We can now start estimating. First, for a given C parametrized by (U, c), Kitaoka proved (see \26\ 
Prop. 1]) that the Kloosterman sum satisfies 

|^(Q,r;C)|« C ?cf +£ (c 2 ,TH) 1 / 2 , 

for any e > 0, where v is the second column of V and the implied constant depends only on e. 
Hence by (|3.1.1ip . we obtain 

£ hk(M, T) « (LM) fe - 3 / 2 cl /2 c 2 - 1+£ (c 2 ,rM) 1 /2| ajfc (p( C ))|. 

D (mod C) 

In order to handle the Bessel integral 3k(P(C)), we will partition M|(Z) in three sets Ci, C 2 , 
S3, according to the relative sizes of the values s\ and s 2 for the corresponding invariants c. These 
can be determined from the size of Tr(P) and |P|; precisely, we let 

d = {C I Tr(P) < 1}, 

C 2 = {C I Tr(P) >max(2|P|,l)}, 

C 3 = {C I 1 < Tr(P) < 2|P|}, 

and we further denote by Gj(c) the subsets of where C is associated with the invariants c = 
(ci,c 2 ,y). The following lemma gives the rough size of these sets, or a weighted version that is 
needed below: 

Lemma 3.4. With notation as above, for any c = (ci,c 2 ,y), we have 
(3.1.15) |Ci(c)| < (ac)-^ 2 - £ , 

(3-1.16) £ |A|^(Ti(AI^))- 8 /^<(^^ ^ C <| 

(3.1.17) |e 3 (c)| « (ac) 1 / 2 ^, 

/or any e > and S > in the second, where the implied constants depend on 5 and e. 

Proof. All these are proved by Kitaoka. Precisely: 

• The bound (13, 1.15ft comes from |26} p. 167], using the fact that in that case we have a <C 1; 

• The bound (13. 1.16ft comes from the arguments of [26, p. 168, 169] (note that in that case 
the summation set is infinite); to be more precise, Kitaoka argues with what amounts to 
taking 

6 = k/2- 7/4, so that fc/2 -3/4=1 + 5, (1 - k)/2 = -5/4 - 5, 

but the only information required (up to [26, p. 169, line 10]) is the sign and the value of 
the sum of the two exponents 

k/2 - 3/4 + (1 - k)/2 = -1/4 = 1 + 5 + (-5/4 - 5) 

(this is used in [261 P- 168, line -12]). Hence Kitaoka's argument applies for 5 > 0. 

• The bound (13.1.17ft comes from [26, p. 168], using the fact that in that case we have c> 1. 
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□ 

As shown also by Kitaoka, we have the following crucial localization properties (see \26\ p. 166]): 

(a) If C e Ci, then s\ < 1 and s 2 < 1; 

(b) If C € 62, then s\ < 1 and S2 S> 1, with absolute implied constant; 

(c) If C £ 63, then si 3> 1 and S2 3> 1, with absolute implied constants. 

Now, by breaking up the sum over C in R2 according to the three subsets Sj, we can write 

I i?2 1 <C -R21 + -^22 + R23 
where, for i = 1, 2, 3, and any fixed e > 0, we have 

i? 2l « (LM) fe ~l c i c 2 1+e ( c 2, rfu]) 1 / 2 ^^), 
cev 

for any e > with 

^(c)= 2 l^(C))|, 

cee,(c) 

the implied constant depending only on d and e. 

Accordingly, we study each of i?2i, i?22, -R23 separately. 

— Estimation of R21. Since we have s\ < 1,S2 < 1, we use (|3.1.3p : using the superexponential 
growth of the Gamma function, we obtain easily 

\UP(c))\ «s 

for C £ Ci and any fixed 5 > 0. On the other hand, by (|3.1.15p . we have 

|Ci(c)| < (ac)- 1 ^ « (si S2 )-^ 2e , 



(3.1.18) * 2 i(c) « « {LM)^^%- 



for any e > 0, and taking it small enough we obtain 

2 

for any fixed 5 > 0. For fixed c\, C2 first, we have 

y6SL(2,z)/ro( C2 /ci) v 
from which one deduces easily 

£ cFc 2 - 1+£ ( C 2,TM)V^ 21 ( c) « ^ (ci ^( C2 / Cl)i M 2 )V 2 

veSL(2,z)/r°(c 2 / Cl ) 

for any 5 > 0, possibly different than before (using (|3.1.13p and [26j Prop. 2] to handle the gcd; 
the exponent of c\ was worsened by 1/2 to facilitate the use of this lemma). 

Writing C2 = nc\, with n ^ 1, we can finally sum over c\ and n; the resulting series converge for 
5 > and we obtain 

R21 « 2- k (LM) k -'i +l+s Y, cfn-^+'i^LM 2 ) 1 / 2 , 

ci,n^l 

and therefore by taking, e.g., 5 = 2e (and changing notation), we derive 

(3.1.19) R 2 i < (LM) k -l(LM) l+£ k- E 

for any e > and E ^ 1, where the implied constant depends on d, E and e. 
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— Estimation of i? 22 . We treat the R22 term next. Using (|3. 1 .3|) for the Bessel function 
involving s\ and (|3.1.5p for the one involving s 2 , and using the fact that 

Tr(A[U]) x Tr(P) = s\ + s 2 2 X sj 

for this term, it is easy to check that 

|9fc(jP(C))l T(k-3/2) ^" i{TrA[U])1 ^- 

If we write 
for any fixed 5 > 0, and observe that 



\A\ ac 



Tr(A[U}) ~ s\ * Sl 1; 
it follows using the super-exponential growth of the Gamma function that 

ft 22 (c)«2- fc \A\ 1+5 (TrA[U})- r i- s 
c*ee 2 (c) 

for any fixed 5 > 0. By (|3.1.16p . we have 

[(ac)4 +e if ac > 1, 

for any e > and 8 > 0. 

We take e = 5/2 and using (|3. 1. 14[) . we deduce that 

V cic 2 >d 1 /2iM/4 V 2 7 VeSL(2,Z)/r°(c 2 /ci) 



+ E fe^" E ck-c*.™* 

cic 2 <d 1 / 2 iM/4 \ / FeSL(2,Z)/rO(c 2 / Cl ) 



for any 5 > 0. In the second sum, we can write trivially 

( d^LM \ l+& < / dV 2 LM \ 1+<5 <<; / LM N 1+5 
y 2cic 2 y y 2cic 2 J \cic 2 

so we end up with 

%«^EO W E ck-WM)- 

ci|c 2 V 7 VeSL(2,Z)/r°(c2/ei) 

and now, using the same type of arguments leading from (|3.1.18|) to (|3.1.19p . we see that 

(3.1.20) R 22 « (LM) k - 3 / 2 (LM) 1+£ k- E 

for any e > and i? ^ 1, the implied constant dependind on d, E and e. 

— Estimation of i? 2 3. Recall that we have 1 <C s\ < s 2 for C £ S3. We estimate the Bessel 
integral using 



\3k(P)\^( + ) J fc _3/2(47TSi sin0) J fc _3 /2 (47r52 sin^) sin 6 



d0 
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where 

Mi = {6 G [0,7r/2] | 47rsisin6l < 1}, 
M 2 = {0 G [0,vr/2] | 1 < 47rsisin0, 1 < 47rs 2 sin0}. 

In the first we use (|3. 1 .3|) and the super-exponential growth of the Gamma function to write 

4-3/2 (4ttsi sin 0) < 2~ k s{, 4_ 3/2 (47rs 2 sin0) 

for any 5 > 0, and in the second we use the estimate (|3.1.4p to get 

Jfc_3/2(47rsisin0)j fc _3/ 2 (4ffsasin0) < *T 2/3 , 

so that 

|0 fe (P(C))| « AT 2 / 3 + 2~ fc ( SlS2 )' 5 « Ar 2 / 3 ( Sl s 2 ) 5 
for any 5 > 0. It follows that 

3Wc)«fc- 2/3 ^ (s lS2 ) 5 , 
c*ee 3 (c) 

which, by (|3.1.17p with, e.g., e = 5, gives 

^ 23 (c) < fc- 2 / 3 (LM) 1+5 (c 1 c 2 )- 1 - 5 

for any 5 > 0. Then the same argument as that following f|3. 1. 18[) is used to sum over the parameters 
c, and to deduce 

(3.1.21) R 23 < (LM) k - 3 / 2 (LM) 1+£ k- E 

for any e > and E ^ 1, the implied constant depending on d, S and e- 
Summarizing, we have 

(LM) k ~ 3 / 2 A(k; c, c', L, M) = R l + R 2 ^R l + R 21 + # 22 + R 23 , 

and putting together the estimates (|3.1.9p . ()3. 1. 19|) . (|3.1.20p . and ()3.1.21|) . we find that we have 
proved the estimate 

A{k; c, c , L, M) < L 1+e M 3 / 2+£ AT 2 / 3 
for e > 0, which was our goal. □ 

Remark 3.5. For later investigations, it may be worth pointing out that the limitation on the error 
term, as a function on k, arises only from the contributions R\ and (the second part of) R 2 %. All 
other terms decay faster than any polynomial in k as k — > +00. 

3.2. A quasi-orthogonality relation for Siegel modular forms. We now put together the 
results of the previous sections. For every k ^ 1, we fix a Hecke basis §* k of Fix the data (d, A) 
as in Section [L2l and let w^ rfA be as defined there; accordingly we have measures v sop ,k defined for 
every finite set of primes S and weight k ^ 1 using suitable average over F G S£. 
Our main result in this section is: 

Proposition 3.6. Let S be a finite set of primes, and I = (l p ), m = (m p ) be S-tuples of non- 
negative integers. Put 

L = J[p l ", M=Y[p m i>. 

pes pes 

Then we have 

f J] U^ mp {x p )du s ,k = £ < d , A II ujr' mp (a p (F), b p (F)) — ► 5(1; m) 
JXs P es Fes* pes 
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as k —> oo over the even integers, where 
5(1; m) ■- 

More precisely, for any even k we have 



1 if L = M = 1, i.e. all l p and m p are 0, 
otherwise. 



' L l+e M 3/2+e' 



(3.2.1) Yl <d,A II U l p p ' mp (a p {F),b p (F)) = 5(1; m) + O 

Fes* pes 

for any e > 0, where the implied constant depends only on d and e. 

We will first prove a lemma which is easy, but where the distinction between SL(2, Z) and 
GL(2, Z)-equivalence of quadratic forms is important. 



Lemma 3.7. For c, d G CLj, put 

5(c,c') = 

Then we have 



1 if c is GL(2,Z)- equivalent to c' , 
otherwise. 



Y, A( C )A^( C , C 0|AU,( C )| = ?MZ^M) 
c,c'ec\ d A 

where 



d A 



1 z/A 2 = l, 

2 otherwise. 



Proof. Let H C Cl^ be the group of 2-torsion elements. The classes d which are GL(2, Z)-equivalent 
to a given class c are c and c , hence there are either one or two, depending on whether c G H 
or not. Similarly, | Aut(c)| (which is the order of GL(2, Z)-automorphisms of a representative of c) 
equals either 2w(—d) or w(—d), depending on whether c lies in H or not. 
Therefore, we have 

£ A( C )A^( C , C ')|Aut(c)|= £ A(c)\Ant(c)\J2W)S(c,c') 

e,c'eCl d eeCl d d 

= Y I Aut(c)|A(c) 2 + ^ I Aut(c)|A(c)(A(c) + A(c -1 )) 
= «,(-d)£(l + A 2 (c)) 

c 

by writing 2 = 1 + A(c 2 ) = 1 + A 2 (c) when c £ H. The result follows immediately. □ 

Now we come to the proof of Proposition 13.61 
Proof. For brevity, we drop the subscripts d and A here. For F £ S£, by definition of oj£ , we have 



, 3 



s C/ / p (a p (F),6 p (F)) = — p^-^ — [[US p (a p (F),b p (F)). 



Now, we write 



\a(d, A; = a(d, A; F)a(d, A; F) 
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and using (jl.2.3p to express the first term, we get 

|2 



Hd ( p^ 1 II U l p p ' mp (a p (F),b p (F)) = ^ WHF,S c ) n U l p ^(a p (F)MF)). 

^ ' ' pes ^ ' ' ceci d pes 

Now Theorem 12.101 applies to transform this into 



L,M\ 



a(d,A;F)\ 2 n i p ,m p[ fm , fjA . \Cl d \ L2~ k M 2 ~ k a(d, A; F) ^ L 
v ' ' pes \ ' / deC\ d (M) 

_\C\ d \L^M^ ^ Hc)W) W^T)a{F,S L c ; M ) 



CUM)\ ^ (F, F) 

c'eC\ d {M) 



after expanding a(d, A; F) using its definition. We are now reduced to a quantity involving only 
Fourier coefficients. 

We then apply the basic property of the Poincare series (13.1.ip to express these Fourier coefficients 
in terms of inner product with Poincare series: we have 

-fc+3/2 



(F,P kySc ) = 8c k (-) a(F,S c ), 

/ ,\ -fc+3/2 

(F, P ks l M ) = 8c fc (LM)" 2fc + 3 ( | J a(F, S 1 / 1 



for c G CLj, c' G Clrf(M), and multiplying out with the normalizing constants, we get 



M k - 1 L k ~ 3 / 2 d A (d/4) k -l ^ .^( F i P k,s c }(F,P k>s L,M) 



c'eci d (M) 

for every F G S£. 

We now sum over i 7 and exchange the summation to average over F first. Since {F/ll^ll} is an 
orthonormal basis of the vector space § k , we have 



\F n 

Now, according to (|3.1.ip again, we have 



/dL 2 M 2 \-k+l 
(Pk,S c ,P kiS L,M) = 8c k{-^—) a(k;c,c',L,M) 

where a(k; c,d ,L, M) denotes, as before, the S^' M -th Fourier coefficient of the Poincare series Pk,s c 
Applying this and the formal definition 

a(k;c,c',L,M) = | Aut(c)|<5(c, c'; L; M) + L k ~ 3 / 2 M k ~ 2 A(k; c, c', L, M), 

as in Proposition 13. H we obtain first, using Lemma 13.71 that 

£ u>l J] U l p v,rrip (a p (F) ,b p (F)) = 5(1; m)+ 
Fes* pes 

^ +3/2 ^ t+ % M - dW - d) |CU(M| E Mc)W)M^.L,M), 

c'ec\ d (M) 
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and then Proposition 13.11 and Proposition 13.31 lead immediately to the desired result. 



□ 



Remark 3.8. In the case of cusp forms on SL(2, Z) and its congruence subgroups, one can write the 
Petersson formula in a way which is suitable for further transformations (with "off-diagonal terms" 
involving Kloosterman sums), as first investigated by Duke, Friedlander and Iwaniec. These are of 
crucial importance in, e.g., the extension of the range of test functions for low- lying zeros in [24J). 
In our case, the complexity of the analogue expansion (which is only implicit in Kitaoka's work) 
for Siegel cusp forms makes this a rather doubtful prospect, at least at the moment. 



4. Local equidistribution 

To pass from Proposition 13.61 to a local equidistribution result, we must understand how the test 
functions considered there relate to the space of all continuous functions on Y5. This is the purpose 
of this section. 

4.1. Symmetric functions and polynomials. We first observe explicitly that the Laurent poly- 
nomials 

U l p > m (a,b) G C[a,6,o- 1 ,6- 1 ] 
of Theorem 12.51 are invariant under the transformations 

(a, b) i-)- (b, a), (a, b) H> (a -1 , b), (a, b) H> (a, 6 _1 ), 

which means that they can be interpreted as functions (also denoted U p m ) on the space Y p or 
on the set X p of unramified principal series of G(Q P ). We first state a simple consequence of 
Proposition 12.71 

Corollary 4.1. Let S be a fixed finite set of primes, and let Yg be as before. The linear span of 
the functions 

(a p ,bp) p eS •->• JJ U l p p,mp (a p ,b p ), 
pes 

where (l p ), (m p ) run over non-negative integers indexed by S, is dense in the space C(Y$) of 
continuous functions on Y$. 

Proof. By the Stone- Weierstrass Theorem, this follows immediately from Proposition 12. 7\ using the 
product structure to go from a single prime to a finite set of primes. □ 

The point of this, in comparison with Proposition 13.61 is of course the following fact: 

Proposition 4.2. Let S be any fixed finite set of primes, and let /is be the associated Plancherel 
measure on Ys, defined in the introduction. We have 



f TT TTh^v t u \ a J 1 iflp = m p = for all p£ S, 
/ U^p p (a p ,b p )dfi s = i _ • 

J y s I U otherwise, 



pes 

for all non-negative integers (l p ), (m p ) indexed by primes in S. 

Proof. Since we work with product measures and product functions, it is enough to prove this for 
the case n = 1. But that follows directly from |17l equation (8)]. □ 

Remark 4.3. This fact can also be proved by direct contour integration via Cauchy's formula using 
the generating function description for U p m (a,b) (given by Theorem 

It is now a simple matter to conclude the proof of Theorem 11.61 

29 



Proof of Theorem \1.6l Fix a finite set of primes S. Using the Weyl equidistribution criterion, in 
order to prove that vs,k converges weakly to us as k — > +00 over even integers, it suffices to show 
that 

(4.1.1) lim / ip((x p ))dv Sk = \ tp(x)d\x s {x) 

j Ys Jy s 

for all functions (p taken from a set of continuous functions whose linear combinations span C(Ys). 
By Corollary 14. 1\ the functions 

¥>((%>, M) = II Up p,mp (a p ,b p ), 
pes 

where (l p ), (m p ) are non- negative integers indexed by S, form such a set. But for 93 of this type, 
the desired limit (|4.1.ip is obtained by combining Proposition 14.21 and Proposition 13.61 

Now, for the proof of the quantitative version f 1 1 . 2 . 6 1) . First of all, we can assume that all 
polynomials (p p are non-constant, i.e., that d p ^ 1 for each p G S, by working with a smaller S 
if necessary (and incorporating the constant functions at a single prime). The polynomials (p p are 
finite linear combinations, say 

ip p (a p ,b p ) = ^ E 0p( l pi m p) u p ,mp ( a i b ) 

of the basis polynomials U p p ' mp , for some e p , f p ^0 with max(e p , f p ) ^ 1. 

Taking the product of these expressions over S, summing over F, and using fj3.2. If) we get 

(4.1.2) Yl «>L<P(M F )M F ))pes) = [ ^(x)dfi s (oo) + k- 2 / 3 R 

where we see that the remainder R can be bounded by 

1^1 «E E L^M^J\\^ p (v p (L\v p (M))\ 

L\L V M\M V peS 

for any e > 0, where the implied constant depends only on e and 

L V =n P ep , M ip =i[pf*. 

pes pes 
The coefficients in the expansion are obtained as inner products 

<p p (!,m) = &p,U l p > m ) 

in L 2 (Y p ,diJ, p ), by orthogonality of the polynomials U p m . Since the underlying measure \x p is a 
probability measure supported on the tempered subset X C Y p , those coefficients may be bounded 
by 

\<p p (v p (L),v p (M))\ ^ ||^t/p p ' mp ||oo ^C( ? ; p (L)+3) 3 (t; p (M) + 3) 3 ||^|| 00 , 
by Lemma 12.81 Therefore, we get the estimate 

\R\ « L^M^rj(LM) J] = L^Mf^LM) |MU 

pes 

where 77 (n) is the multiplicative function such that 

for p prime and v ^ (here we use the fact that max(e p , f p ) ^ 1 for each p). This is a divisor-like 
function, i.e., it satisfies 

7?(n) <C n s 
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for any e > 0, the implied constant depending only on e. Therefore we get 

for any e > 0, where the implied constant depends only on e > 0. 

To derive (|1.2.6p . we observe that, by the second part of Proposition ^, 7\ the linear decomposition 
of (fp holds with 

Gp ~\~ 2fp dp, 

where d p is the total degree of <p p as polynomial in (a + a _1 , b + b^ 1 ). Thus, at the cost of worsening 
the factor involving M, we obtain (jl.2.6p , □ 

Remark 4.4. The proof shows that if we know that the factors (p p are combinations of polynomials 
Up" 1 with I ^ l p , m ^ m p , we have the stronger estimate 



V ooi kV ((a p (F),b p (F)) P es) = [ <p(x)dns(x) + ofk^L^M 3 / 2 ^ 
for any e > 0, where 



pes pes 
5. Applications 

We now gather some applications of the local equidistribution theorem. To emphasize the general 
principles involved, and their expected applicability to the most general "families" of L-functions, 
we denote 

B k (a(F)) = —L^ £ ^a(F) 

2^ u k Fes* 
Fes* 

for k ^ 2 even and for any complex numbers (a(F)). This is the averaging operator for a probability 
measure depending on k, and we know from the previous results that 

E fc (a(F)) ~ ]T u£a(F) 

Fes* 

as k — > +oo over even integers. We denote by P&(») the associated probability. We also recall 
that §2k 1S the se t °f CUS P forms which are not Saito-Kurokawa forms, and we denote by SS^ the 
complementary set of Saito-Kurokawa lifts. 

5.1. Direct applications. We start with direct consequences of the local equidistribution. The 
first is partly superseded by the proof of the generalized Ramanujan conjecture in our case [50] . 
but it may still be taken as an indication that the special Saito-Kurokawa modular forms which 
fail to satisfy it are "few", even when counted with our special weights. 

Proposition 5.1. (1) Fix a prime p. Then "most" F € § k satisfy the generalized Ramanujan 
conjecture at p, in the sense that we have 

f 1 2k('Kp(F) is not tempered) = P2k{^p(F) ^ X C X p ) — > 

as k — )• +oo. 

(2) Let a(F) be any bounded function defined for F which are Saito-Kurokawa lifts, i.e., F G §|. 
Then we have 

lim B 2k (a(F)l , ) = 0. 
Proof. Since the limiting measure \i p is supported on X C X p , this is immediate. □ 
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In particular, it follows that the measure v s . defined as vs,ki Du t with F restricted to § k , also 
converge weakly, as k — > +00, to us- We will denote by ~E k (a(F)) the average 

Ei(a(F)) =E k (a(F)l {Fegi} ), 

so that the previous result means that this is still, asymptotically, a probability average. 
The next result has the feel of a "Strong Approximation" theorem: 

Proposition 5.2 ("Strong approximation"). (1) Let Aut denote the set of all cuspidal automorphic 
representations on GSp(4, A). Then for any finite set of primes S, the local components tts for 
those it € Aut unramified at the primes in S form a dense subset of X s . 

(2) Fix a finite set of prime S, and let (e p ) pe s be signs ±1. There exist infinitely many Siegel 
cusp forms F of level 1 which are Hecke-eigenforms such that the Hecke eigenvalues at all p G S 
have sign e p . 

Proof. (1) The support of the limiting measure /is is X s , hence the result is again immediate 
(with the much more precise information that denseness holds already for it associated to Siegel 
cusp forms of full level). 

(2) This sample application follows from the fact that the Hecke eigenvalue at a prime p is 
a p + a~ l + b p + b~ 1 for F 6 S£, and it is clear from the formulas for the Haar measure \i (|1.2.4p 
and for the density function A p (jl.2.5j) that 

fj. p ({(a p , b p ) peS I the sign of a p + a~ 1 + b p + bp 1 is e p }) > 

for any prime p and hence 

Ms({(%>i bp) p( zs I the sign of a p + a" 1 + b p + b' 1 is e p for p e S}) > 

(and hence, for k large enough, at least one F £ §>* k satisfies those local conditions). □ 

Remark 5.3. In fact, our local equidistribution shows much more. For instance, for any non- 
negligible subset of T of X s and any fundamental discriminant —d, our result shows that one can 
find infinitely many Siegel modular forms whose local components at S lie in T and the sum of 
whose Fourier coefficients of discriminant — d is non-zero. Since such sums of Fourier coefficients 
of Siegel modular forms are conjecturally related to central critical L values of the twisted forms 
- see Section 15.41 — this can be interpreted as a (conditional) result on the plentitude of Siegel 
modular forms with prescribed local behavior and non-vanishing central critical values. 

5.2. Averaging L-functions. The spin L-function associated to F G £>* k is defined, in terms of 
the Satake parameters (a p ,b p ), by the Euler product 

l(f, s) = J] (1 - a pP - s )-\i - bp P - s )-\i - ^ V s ) _1 (i - K 'p-T 1 ; 

p 

it is explained in [2] that this is a Langlands L-function associated with the Spin representa- 
tion Spin(5,C) — > GL(4, C). This Spin group is the Langlands dual group of Sp(4), and since 
Spin(5) ~ Sp(4), this is natural in view of the parametrization of the local representations in terms 
of semisimple conjugacy classes of USp(4, C) which is described in the introduction. From the 
point of view of Sp(4), this is the Langlands L-function corresponding to the natural representa- 
tion GSp(4) C GL(4). 

The idea is that in the region of absolute convergence (which is Re(s) > 1, for F not a Saito- 
Kurokawa lift), the average of such a product is the average of asymptotically independent random 
variables, and hence will be the product of the averages for the limiting distributions at each p. 
Saito-Kurokawa lifts, being asymptotically negligible, do not cause much trouble in that case. 

To go to the details, we first recall that, from work of Andrianov [lj, it is known that L(F,s) 
has the basic standard analytic properties expected from an L-function; it is self-dual with root 
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number (—1) and unramified at finite places; it extends to a meromorphic function of s, and the 
completed L-function 

A(L, s) = (2tt)- 2s F(s + l/2)r(s + k — 3/2)L(L, s) 

satisfies 

A(F,s) = (-l) k A(F,l-s). 
Furthermore, if k is odd, the L-function is entire; otherwise, it may have poles at s = 3/2 and 
s = —1/2, and this happens precisely when F is a Saito-Kurokawa lift (i.e., if F 6 sjj. using the 
notation of Theorem II .ip . In that case, the L-function is given by 

(5.2.1) L(F, s) = C(s - l/2)C(a + 1/2)L(/, s) 

for some classical cusp form / of weight 2k— 2 on SL(2, Z) (the L-function of which is also normalized 
so that the critical line is Re(s) = 1/2). Note that although, in general, there are other automorphic 
forms on GSp(4) where the L-function have poles, the Saito-Kurokawa lifts are the only ones which 
are holomorphic with level 1 (we refer to Piatetski-Shapiro's paper [36] for more details). 

Remark 5.4. For completeness, even though we do not need it here, let us recall the corresponding 
results for the other "standard" L-function, which is the degree-five L-function coming from the 
projection 

pr : Spin(5,C) -»• SO(5,C) C GL(5,C). 
This L-function has the form 
(5.2.2) 

L(F, pr, S ) = n ((1 " " a p b pP - s )(l - a^p-^l - a p %p- s ){l - (opfep)" 1 ?-')" 1 - 

p 

From work of Mizumoto [M], it is known that L(F, pr, s) has the basic standard analytic proper- 
ties expected from an L-function; it is self-dual with root number 1 and unramified at finite places; 
it extends to an entire function of s, and the completed L-function 

A(F, pr, s) = 2- 2s ir- 5s/2 T (^y^) T(s + k- l)T(s + k- 2)L(F, s) 

satisfies 

A(F,pr,s) = A(F,pr,l-s). 

Let us now return to the spin L-function. To average the L-function, we express it in additive 
terms. For this purpose, we denote by 

7T p : § k — > X p 

the map F \— > ir p (F), the local component of the automorphic representation np associated with F 
as described earlier, which we identify with (a p (F) , b p (F)) E Y p . 

Expanding the Euler factors in powers of p~ s and then expanding the product into Dirichlet 
series, we find the expression 

L(F,s) = J2KF,n)n- s 

in the region of absolute convergence, where 

X(F,n) = HH Mn) (7r p (F)), for n = \{p^ n \ 

p\n p\n 

in terms of functions H m , m ^ 0, on Y p given by the symmetric functions 

H m (a,b)= a kl ~ ks b k2 - ki 

ki+k 2 +k 3 +k4=m 
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(note that H m is independent of p, though that is not crucial in what follows, and that it is 
well-defined on X p since it is invariant under the Weyl group). 

If Re(s) > 1, the series L(F, s) converge absolutely at s for F G § 2k , and we have 



E^(L(F, S )) = ^E^(A(F,n))n 
Fix n first, and factor it as before 



n = Y\v Vv 

p\n 



By our local equidistribution theorem applied to u s 2fe , we have 



Et> 
2k 



(\(F,n))=El k (j[H Vp (Tr p (F)f)^llj H Vp (x)d^{x) 

p n p\n 



as k — > +00. Therefore, by the dominated convergence theorem, we have 

t (x)d(ip(x))n~ 



ri^l n^l p\n Jx p 



since, using the formula 



~(m + l)(m + 2)(m + 3) 



for the number of monomials of degree m in 4 variables (the number of terms in H m ) we have 

E^ k (X(F,n))n- s \^n-ll(v p + 3f 

p\n 

for all n ^ 1 and fc, which defines an absolutely convergent series for a = Re(s) > 1. (We use here 
the generalized Ramanujan conjecture, proved in this case by Weissauer [50].) 
Now we refold back the limiting expression as an Euler product: 



e (n / H VV (x)dn P {x))n- s = nEf- ,a / W^Mp 

= II / L p(x,s)dfj, p (x), 



where L p (x, s) is the local L-factor of a local representation x = (a, 6) E Xp defined in the lemma 
above (the Euler expansion is justified again by the fact that the series on the left is absolutely 
convergent, as we checked in the lemma). Thus we have proved 

(5.2.3) lim E 2k (L(F,s)) = TT / L p (x,s)dfj, p (x). 

k— >+oo IX 

p J s>-p 

Now assume Re(s) > 1 and s 7^ 3/2. Then all spin L-functions of Saito-Kurokawa lifts are well- 
defined at s, and we therefore also want to have average formulas involving them. If Re(s) > 3/2, 
this is immediate by the previous argument. Otherwise, we have 

^ k L(F,s) = ((s-l/2)((s + l/2) ^kHfF,s), 
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where fp is a classical modular form from which F arises. The L-function L(fp, s) is now absolutely 
convergent, and its values are bounded for all Saito-Kurokawa lifts (by the generalized Ramanujan 
conjecture, for instance). Thus we have 

lim ^2 k L{F,s) = Q 

fc->-+oo ' — ' 

by Proposition 15.11 (2), and this combined with (15.2.3)) gives the result 

lim E 2 k(L(F, s)) = TT / L p (x,s)dfJ v (x). 

k^+co Y Jx p 

At this point, it is clear how to extend this to other Langlands L-functions. Indeed, let 

p : GSp(4,C) -> GL(r,C) 
be an algebraic representation. The Langlands L-function is defined by 

L(F, p, s) = L(tt f , p, s) = JJ det(l - p(x p (F))p~ s )~ 1 

p 

where 

x p (F^ — Xpi^cip, bp) — 

V ' K l J 

is the semisimple conjugacy class of GSp(4, C) associated with ir p (F). We can expand 
det(l — p(x p (F))p~ s y 1 = H (l-a j (a p ,bp)p- s )- 1 

for some polynomial functions ctj on X p , and then we can repeat the same argument to derive fjl . 1 .3j) 
with 

lim E 2fe (L(F,p,s)) = TT / det(l - p{x p (a, b) )p' s y 1 dfip(a, b), 
fc^+oo p Jx 

when s is in the region of common absolute convergence of all F. 

Finally, to get the precise expression in Theorem 11.11 for the spin L-function, we note that the 
special case (12.3.7P of Sugano's formula (Theorem I2.5f) . with Y = p~ s , gives the explicit decompo- 
sition 

l p (x,s) = (i-\ p p-^-s+ (^)p^y\j2 ul A a ^)p^ 

for any prime p. Applying Proposition 14.2) we get therefore the simple expression 
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L p (x, s)dfi p (x) = 

i - v- 1/2 - a + (f ) p- 1 - 23 

and (using the definition of \ p ) we recognize that this is 

L(A,s + l/2), 

where A is the class group character of % = Q(V — d) defining our fixed Bessel models. When d = 4 
and A is trivial, this is ((s + l/2)L(x4, s + 1/2), which is the formula (jl.l.2|) . 

Remark 5.5. In fact, this second argument for the spin L-function can be used to bypass the first 
one (which therefore requires only that we work with the family of functions U p °(a,b)). 
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Remark 5.6. Note that although Theorem 11.11 was stated in the introduction only for averages with 
respect to the weight lu£ := lu£ d A in the special case d = 4, A = 1, our proof actually works for 
general d and A. 

The proof also gives the following fact concerning the limit averages: 
Lemma 5.7. For p prime, let pip be the limiting measure in the local equidistribution result and let 

L p (x, s) = J] (1 - ap-THl ~ bp-THl ~ a" VT'U ~ b^p'T 1 
v 

be the local L- function for x £ X„. Then the Euler product 



Y[ / L p (x,s)dfip(x) 



is absolutely convergent for Re(s) > 1/2. 

Proof. According to what we have said, we have 

Lp(x,s)dp p (x) 



1 _ App -l/2-» + | L-l-2, 



and the result is then obvious. □ 

5.3. Weights and averages over Saito-Kurokawa lifts. In this section, we will explicitly com- 
pute oj^ := when F is a Saito-Kurokawa lift. This will lead to a stronger version of the 
second part of Proposition 15.11 This simple fact is included because it may be helpful for further 
investigations. 

Let k > 2 be even and let ^2k-2 denote the Hecke basis of the space of holomorphic cusp forms 
on GL(2) of weight 2k — 2 and full level. Let F E 2>* k be the (unique) Saito-Kurokawa lift of 
/ G ^2fc-2> so t na * the spinor L- function is given by (15.2.1j) . As usual, we let 

F(Z) = ^a(F,T)e(Tr(TZ)) 

T>0 

be the Fourier expansion of F. It is well-known (see [13j for instance) that a(F,T) then depends 
only on the determinant of T. In particular, it follows that oj£ dA = whenever A / 1. So, we 
assume that A = 1 and shorten 1 to u£ . 
Let 

/0) = ^c{n)e{nz) 

n>0 

be the cusp form of half-integer weight k — \ on To (4) that is associated to / via the Shimura 
correspondence. Then, by [131 Th. 6.2, equation (6)], we have 

(5.3.1) a(T) = c(d) 

for any positive-definite semi-integral matrix T of determinant d/4. On the other hand, by [4], we 
have 

whenever c{d) is non-zero. 

Using (I5.3.ip . (|5.3.2f) and the definition of oojl , it follows that 

(A8iT) 2 h(-d) F(2k - 3) L(f x Xd , k) 



w {-d){k -l)(k- 2) (4tt)»-3(/, /) L(f, 1) 
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Now we consider the average 

^ T(2fc-3) L(fxxd^) 
J-f (47r)^-3(/ )/ ) L(/,l) • 

J 2k — 2 

It is very likely that one can prove an asymptotic formula for this quantity as k — > +00 (possibly 
using the methods of Ramakrishnan and Rogawski in [40]). However, to deal with it quickly, we 
observe first that L(f x Xd, \) arid L(/, 1) are both non-negative (e.g., because L(f, s) is real-valued, 
has no zero with Re(s) > 1 and tends to 1 as s — > +00, and the ratio is non-negative by the above). 
Then, using the fact that L(f, s) has no Siegel zeros (a result of Hoffstein and Ramakrishnan |20j), 
one gets in the usual way a lower bound 

L(/,l)» J-, 
log k 

and therefore 

^ rgfc-3) UJy. Xd ,±) „ r(2fc-3) 1 

,J- (4.)— (/,/) L(fJ) <<{]0gk) f J? (4^(/,/) i(/Xto 2»' 



Next, from the results of Duke [121 . one gets 



r(2A;-3) 1 

L(/x Xd ,-)«l 



(47T)*-»(/,/> u Aa '2 

for k ^ 2, where the implied constant depends on d. The following Proposition, which strengthens 
Proposition 15.11 (2), is then an immediate consequence: 

Proposition 5.8. Suppose a(F) is a complex valued function defined for Saito-Kurokawa lifts and 
satisfying for some 5 > the inequality 

a{F) < k 2 - 5 

for F e §\ k . Then 

km E fc (a(F)l B ) = 0. 

Using weak bounds, like the convexity bound, this applies for instance to a(F) = L(F, 1/2 + it) 
for fixed t 7^ 0. 

5.4. Low-lying zeros, Katz-Sarnak symmetry type and Bocherer's conjecture. The de- 
termination of the distribution of low-lying zeros of the spin L- functions (assuming the Generalized 
Riemann Hypothesis) for restricted test functions is not difficult once a quantitative equidistribu- 
tion statement is known. Conjecturally, the answer indicates which "symmetry type" (in the sense 
of Katz-Sarnak) arises for the family. However, we will see that the answer in our case is surprising, 
and gives some global evidence towards a well known conjecture of Bocherer. 

We now prove Theorem ll.21 This type of computation is quite standard by now, and is known to 
succeed as soon as "approximate orthogonality" has been proved with a power saving with respect 
to the analytic conductor, which is the case thanks to our quantitative equidistribution theorem 
(precisely, from (|3.2.ip ). We may therefore be brief, as far as technical details are concerned (we 
refer to, e.g., [TT], where families derived from classical GL(2) cusp forms are treated with respect 
to the weight). However, since the main term arising from this computation has some meaning, we 
must justify it carefully. 

As before, note that although Theorem 11.21 was stated in the introduction only for averages with 
respect to the weight uj^ := w£^ A in the special case d = 4, A = 1, we actually prove it for any d 
(we stick to A = 1). 
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Proof of Theorem ] ITU Throughout this proof, denotes v For given F, we write 

T 



the logarithmic derivative of the spinor L-function which is supported on powers of primes and 
where A(n) is the von Mangoldt function and, for n = p m , m ^ 1, we have 

c(F,p m ) = a™ + a~ m + 6™ + b p m = Tr(n p (F) m ), 

where n p (F) is interpreted as a conjugacy class in USp(4, C). 

We can apply the following form of the "explicit formula" (see, e.g., |23t Th. 5.12]) to relate 
sums over zeros to sums over primes involving those coefficients: denoting 

7(s) = (2vr)- 2s r(s + l/2)r(s + k- 3/2) 

the common gamma factor for all L(F,s), for any test-function ip which is even and of Schwartz 
class on R, we have 



p 

We apply this to 



where (p is an even Schwartz function with Fourier transform supported in [—a, a]. After treating 
the gamma factor using the formula 

r' r" 

— (k - 1 + it) + — (k - 2 - it) = 2 log k + 0(t 2 k- 2 ), 

which follows from Stirling's formula (see, e.g., |114 §3.1.1, §3.1.2] for precise details of these com- 
putations) and spelling out the von Mangoldt function, we obtain 

(5.4.1) D ,(F) = / RvW<fa - ^££^ c(W ( mj ||_) +0((logk) 
Averaging over F leads to 

B k (D v (F)) = 0(0) - ^ E E^ E ^ c ^^))K m ioiP)) +0((1 ° gfcrl) - 

As usual, easy estimates give 

lim E E ^MF,P m ))0(rn^) = 

k^+oo log(/c 2 ) ^ ^ p m / 2 KX y " »ry Wfc2W 

(the series over primes being convergent even without the compactly-supported test function). 
In the term m = 1, we have 

E fc (c(F,p)) = E k (a p + b p + ap 1 + fc" 1 ) 

(5.4.2) = B k (U^MF))) + hL = hjL + otf+>hr*l*) 

by (|2.3.6p and (|3.2,ip . and hence the contribution of m = 1, which is given by 

logo / logp \ 



io^py E ^ Ek ^ p ^H]^)) 
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is equal to 
2 



lo# 2 )^ P r Vlog(A; 2 ); Vfc 2 /3logfc p/J 



-2/3+£^ 



for any £ > 0, where 



n/r ( v ^ \- Aplogp / logp \ 



logy \ 1 dy 

2/ [1 ,^Hio^)k#)7 + o(1) 



= 2 / <£(x)dx = </?(0) + o(l), 

./[0,+oo[ 

(since 99 is even), by summation by parts using the Prime Number Theorem and the fact that 
X p = 2 or for primes with asymptotic density 1/2 each, so the average value is 1 (see, e.g., 
Lemma 2.7]). 

Now we consider the term m = 2. Although we could appeal to the general estimate (|1.2.6j) . we 
will use an explicit decomposition and (|3.2.ip . First, using (|2.3.7|) . we have 

(5.4.3) U 2 /(v p {F)) = 1 - ^U^(n p (F)) + c(F,p 2 ) + r(vr p (F)) + O&T 1 ) 

where the implied constant is absolute and 

r(a, b) = 1 + ab + ab' 1 + a~ l b + (aft) -1 
as in Theorem 12.51 By (|2.3.8p . we have 

r(vr p (F))(l + a p ) = U°'\n p (F)) + p a(n p (F)) + 0(p~ l ) 

where the quantities a p <C p^ 1 , j3 p <C p~ x l 2 do not depend on F, and the implied constants are 
absolute. Averaging (with E^(-)) from this last formula and using (I5.4.2p . we find 

E fc (rK(F)))«p- 1 +p 3 / 2+£ A ; - 2 / 3 , 

and from (|5.4.3p . we therefore derive 

B k (c(F,p 2 )) = -l + Oip-'+p^k-V 3 ) 

for any e > 0. Consequently, we see that the term m = 2, after averaging, is given by 

2 \ - logp ^ , lT? 2 ^-/ logP \ _ 2 s^ ^ZP - ( ^EP \ , 
log(P)^ p ^k{c^,p ))<P[\ g ( k2) ) ~ log(P)^ p ^^log(P)J + 

^(TO^ pl+£ + (l0g ' rl ) 
= -N k ((p) + 0((log A;)- 1 + k 2a ~ 2 / 3+£ ), 

where 

Nk{ip) = \ —< 2 ^m) =— + o(1) ' 

by computations similar to that of M((p) before. 

We notice that the contribution of the main terms for m = 1 and m = 2 together are 
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0(0) - ^7T^ = / <p(x)d(TSp(x). 



which gives a main term 



Moreover, the error terms in both are negligible as long as 5a/2 — 2/3 < 0, i.e., a < 4/15. Under 
this condition, we obtain therefore 

E k {D v (F))= [ <p(x)dx + o(l), 

as k — > +oo, which is the desired conclusion. □ 

Remark 5.9. It is interesting to note that, to understand the logarithmic derivative of the spin 
L-function L(F, s), one needs to involve the average of the quantity 

t(tt p (F)) = 1 + a p b p + apbp 1 + a~ 1 b p + (a p 6 p ) _1 , 

which is the coefficient of p~ s in the projection L-function L(L, pr, s) (see f|5.2.2[) ) . This illustrates 
the fact that, in the study of automorphic forms on groups of rank r ^ 2, all Langlands L- 
functions (or, at least, those associated with the r fundamental representations of the dual group) 
are intimately linked, and must be considered together. 

We now comment on the relation of Theorem 11.21 with Bdcherer's conjecture. In the literature, 
a density of low-lying zeros given by the measure dasp (as we proved is the case) is taken as 
an indication of symplectic symmetry type (the basic example being the family of real Dirichlet 
characters). Intuitively, these are families of central L- values of self-dual L-functions with sign of 
functional equation +1 for which the central point of evaluation has no special meaning. However, 
although our families are indeed self-dual, a symplectic symmetry seems very unlikely for our family, 
for at least two reasons: first, 1/2 is a critical point in the sense of Deligne, and second, the forms 
of odd weight have functional equations with sign —1. 

There is a natural explanation for the discrepancy: the Fourier coefficient \a(d, 1; F)\ 2 appearing 
in the weight 

k \a(d,l;F)\ 2 

involves probably deeper arithmetic content than one might naively think. Indeed, in [3j, Bocherer 
made the following remarkable conjecture: 

Conjecture (Bocherer 's Conjecture). For any F G Sj^, there exists a constant Cf depending only 
on F such that for all fundamental discriminants —d < one has 

L(F x X d, \) = Cf- d l - 2k w(-d)- 2 ■ \a(d, 1; F)\ 2 , 

where \d denotes the quadratic character associated to the extension Q(V — d). 

Bocherer proved this conjecture for Eisenstein series and Saito-Kurokawa lifts in [3j. Later, he 
and Schulze-Pillot proved an analogue of this conjecture (for Siegel modular forms with level) in 
the case of Yoshida lifts. More recently, works of Furusawa-Shalika [19J , Furusawa-Martin [16J and 
Furusawa-Martin-Shalika [17J have tried to tackle this problem using the relative trace formula. 

Bocherer did not make any speculation about the value of the quantity Cf- However recent 
works such as [16] give some inkling of what to expect. 

We now show that a certain assumption on Cf explains our result on low- lying zeroes. To be 
more definite, we make the following hypothesis. 
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Hypothesis. For non-Saito-Kurokawa forms F £ S^, we have 

(5-4.4) uf = L(F, l)L{F x Xd , \)L{ Xd , irHF) 

in terms of spinor L-functions, where j(F) > is "well-behaved" , in particular 



7(F) 



has a positive limiting average value as k — ^ +00, and j(F) is asymptotically independent of the 
central special L-values. 

In terms of Fourier coefficients, this hypothesis is equivalent to the following specific variant of 
Bocherer's conjecture: for all F 6 S^fc that is not a Saito-Kurokawa lift, we should have 

(5-4.5) L(F, \)L{F x Xd , I) = ^c 2kl {F)-\d/A)^ k W (-d)-' W £^ ? • 

Remark 5.10. Such a formulation (involving two central values, or in other words a central value for 
the base-change of F to the quadratic field Q(V~ d)) is strongly suggested by [16j (1.4)] and |39j. 
It is also compatible with a conjecture of Prasad and Takloo-Bighash [39], which itself is an 
analogue (for the case of Bessel periods) of the remarkable Ichino-Ikeda conjecture |21j dealing 
with (SO(n), SO(n — 1)) periods. In this context, it is also worth mentioning that the question of 
vanishing of the Bessel period, i.e., the vanishing of a(d, A; F), is closely tied with the Gross-Prasad 
conjecture for (50(5), SO (2)). 

Under our stated hypothesis (|5.4.4p . we consider the crucial average 

E 

for a fixed prime p. Our goal is to show that this allows us to recover naturally the formula (|5.4.2p 
from which the "mock-symplectic" symmetry-type arose in the proof of Theorem 11.21 (the contri- 
bution of p 2 was consistent with the expected orthogonal symmetry type). Thus, assuming (I5,4.4j) . 
we need to compute the average 



E 7(i ? )^(i ? , \)L{F x Xd , \)L{ Xd , iy 1 c(F,p), 



Since c(F,p) = X(F,p) is also the p-th coefficient of the Dirichlet series L(F,s), and since the 
analytic conductor of both L-functions is about k 2 , we see by applying a suitable Approximate 
Functional Equation (and recalling that the sign of the functional equation is 1 for both L-functions) 
that this is roughly 

2L( Xd ,l)- 1 E ^0 E l(F)X(F,m)HF,n)X(F,p) 

(the sums should involve a smooth cutoff). 

Under the (reasonable) assumption that the coefficients of the Dirichlet series are asymptotically 
orthogonal under this average^ one is led to the guess that the terms which contribute are those 
with m = np or n = mp, and thus one should have 

E J(F)L(F, \)L(F x Xd , l)L( Xd , l)^c(F,p) « 1±*M = ± 

F % 2 2 VP VP 



This depends on the hypothesis that 7(F) is innocuous. 
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as k — >■ +00, where the L(l,Xd) has cancelled out with 



\- Xd(m) 



in 



(again with a smooth cutoff). 

But this is exactly what we proved in (|5.4.2|) . and what led to Theorem 11.21 We therefore 
interpret this as a (global, averaged) confirmation of Bocherer's conjecture in the form (I5.4.5p . 



Appendix: comparison with GL(2)-families 

This section is intended to summarize some basic facts about holomorphic Siegel modular forms 
and their adelic counterparts, by comparison with the case of classical modular forms for congruence 
subgroups of SL(2,Z). We also give references for the SL(2)-analogues of the results in this paper. 

• The closest analogue of our family of cusp forms is the set Jit of primitive holomorphic cusp 
forms of weight k for SL(2,Z), with trivial nebentypus, counted according to the weight 
given by 

_r(fc-i) 1 
">/- (47r)fc _i {fjy 

In contrast with §£, this set is the unique Hecke-eigenbasis of the space "Kt of cusp forms 
of weight k and level 1; in our context, the corresponding multiplicity one theorem is not 
known (because the Fourier coefficients are not functions of the Hecke eigenvalues), and so 
the Hecke basis S£ is not necessarily unique in 

Another obvious distinction is the presence of the Fourier coefficient a(F, 1) in (jl.l.ip . 
As we saw, this has crucial arithmetic content. A way to insert this aspect "by hand" into 
the classical case is to consider the twisted weights 

u' f = au f L(f, 1/2), or J f = au f L{f, 1/2)L(/ x X d, 1/2), 

where a > is a constant such that 

as k — > +00. (The existence of the limit that makes this normalization possible is essentially 
already in Duke's paper |12} Prop. 2], where the limit is with respect to the level.) 

• The local equidistribution theorem for "K* k , as k — > +00, is the following: for any prime 
p, the map sending / to the p-component of the associated automorphic representation of 
GL(2, A) can be identified with the map 

P ^ A/(p) € R 

k — 1 

where p~^~ \f(p) is the p-th Hecke eigenvalue, or equivalently the p-th Fourier coefficient. 
By Hecke's bound, we have \f(p) € [—2^/p,2y/p], and the associated representation of 
GL(2, Q p ) is the unramified principal series obtained from the unramified characters a, /3 
of Qp such that 

a(p)P(p) = 1, a(p) + p{p) = X f (p). 

Then, for any finite set of primes S, the measures ms,k defined as the sum of Dirac 
measures at A/ (p) for p G J{£ converge weakly to the measure 

ns = JJ fJ-ST 
pes 
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where fisT is the Sato- Tate measure, supported on [—2,2], given there by 

2 

— yl - x 2 /4dx. 

TT 

The above fact is quite easy to prove. First, the Hecke relations describe Aj(L) in terms of 
the factorization of L ^ 1, namely 

\ f (L) = l[u lp (\ f (p)), 

p\L 

where U[ is the l-th Chebychev polynomial defined by 

n , sm((Z + l)0) 
!7 z (2cos0) = U . - 1 
smt) 

These form a basis of polynomials in one variable, and hence span a dense subset of 
C(HVP, 2-^), with 



Ut(x)d» ST (x) = 5(1,1). 
The second ingredient is the Petersson formula; indeed, for any l p ^ 0, let 

L = ]Jp lp , 

pes 

and then we have 

E ^/IW A /(p))= E w / A /( L ) 

f&* k pes feK* k 

= 8(L, 1) - 2vrr fc C ~ 1S ( L > 1; c ) J fc-i (^?) — »• 5(L, 1) 

as k — > +oo, where S(L, l;c) denotes the standard Kloosterman sum. This gives the local 
equidistribution statement. Note that, in contrast with our results, the limiting measure at 
p is independent of p in this case. 

The Hecke relations are analogues of Sugano's formula (Theorem [23]); the Petersson formula 
and the related orthogonality are the analogues of Proposition 13.61 On the other hand, 
the necessary work to go from Fourier coefficients (controlled by Poincare series) to Hecke 
eigenvalues is completely absent from the classical case. 

Analogues of the direct applications of Section [57Ll were proved first, essentially, by Brugge- 
man [5] (analogue of Proposition 15.11 for Maass forms, where the Ramanuj an- Petersson 
conjecture is not yet known); analogues of Proposition 15.21 are due to Sarnak [33] (Maass 
forms) and Serre (holomorphic forms), though both counted the cusp forms with the natu- 
ral weight 1, and used the trace formula instead of the Petersson formula (correspondingly, 
their limiting distributions was different: at p they obtained the Plancherel measure for the 
unramified principal series of GL(2,Q P ) with trivial central character). 
Computations tantamount to working with the twisted weight lv'j are also classical (in 
particular, the computation of 

E w / A /( m ) 

/ 

for a fixed m is a special case of the first moment computation in |30j . in the case where the 
level goes to infinity and the Rankin-Selberg convolution is the weight 1 theta series with 
L-function Q(s)L(s, xa))- 
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• The analogue of Bocherer's conjecture for "K* k is the famous formula of Waldspurger [49] 
which relates the value L(f x Xdi 5) f° r / ^ t° the squares of Fourier coefficients of 
modular forms of half-integral weight. However, these special values do not appear in the 
standard weights used for averaging L-functions. However, a weighted averaged version of 
Waldspurger's formula was proved by Iwaniec |22| using identities for Kloosterman sums, 
and this may be considered as somewhat analogue to our Theorem II .li 

• The Saito-Kurokawa forms have no analogue in JCt. Indeed, all cusp forms of GL(2) (or 
GL(n)) are expected to satisfy the Ramanujan-Petersson conjecture; for forms in this 
is a theorem of Deligne |10| . On the other hand, Saito-Kurokawa forms do not satisfy the 
generalized Ramanujan conjecture; this is due to the fact that they are CAP representa- 
tions [36] . 
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